NEUTRON  FORM  FACTORS  FROM  ELASTIC 
ELECTRON-OEUTERON  SCATTERING  RATIO  EXPERIMENTS 
AT  VERY  LOW  MOMENTUM  TRANSFERS 


Raymond  William  Berard 


Library 

Nava'  Postgraduate  School 

Monterey,  California  93940 


POSTGRADUATE  SCHOOL 

Monterey,  California 


TH 


NEUTRON  FORM 

FACTORS  FROM  ELASTIC  - 

ELECTRON-DEUTERON  SCATTERING  RATIO  EXPERIMENTS 

AT  VERY  LOW 

MOMENTUM  TRANSFERS 

by 

Raymond 

William  Berard 

and 

Timothy 

Joseph  Traverso 

F.R.  Buskirk 

Thesis  Advisors: 

X.K.  Maruyama 

June    19  7  3 


kppKovtd  Ion  pubtic  kqJLzo&q.)   dUVuhution  uytLuniXfLd . 


T15A906 


Neutron  Form  Factors  from 

Elastic  Electron-Deuteron  Scattering  Ratio  Experiments 

at  Very  Low  Momentum  Transfers 

by 

Raymond  William  Berard 
Ensign,  United  States  Navy 
B,.S.,  United  States  Naval  Academy,  1972 

and 

Timothy  Joseph  Traverso 

Ensign,  United  States  Navy 

B.S.,  United  States  Naval  Academy,  1972 

Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree 

MASTER  OF  SCIENCE  IN  PHYSICS 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
June  1973 


'    ■ 


Libr 

Naval  Postgraduate  School 
Mont, 


ABSTRACT 

Measurements  of  the  ratio  of  the  elastic  electron- 
deuteron  scattering  cross  section  to  the  elastic  electron- 
proton  scattering  cross  section  were  made  for  low  momentum 

transfers.   Augmented  by  the  work  of  Topping  [5],  these  data 

2  -2  -2 

span  a  range  of  q   from  0.05  fm    to  0.50  fm   .   From  these, 

the  electric  form  factor  of  the  neutron,  G„  ,  was  extracted 

hi 
n 

using  three  Lomon-Feshbach  wave  functions  for  the  deuteron, 

2 
(LF1,  LF5,  and  LF15).   Values  of  the  slope  of  G^   versus  q 

hi 
n 

were  measured  to  be  0.0218  ±  0.0022  for  LF1,  0.0220  ±  0.0022 
for  LF5,  and  0.0185  ±  0.0022  for  LF15.   Although  the  data 
seems  to  prefer  model  LF15,  all  were  found  to  be  in  agree- 
ment with  the  thermal  neutron-electron  slope  of  0.019  3  ± 
0.0004.   The  root  mean  square  structure  radius  of  the 
deuteron  was  also  extracted,  using  a  method  devised  by 
Schumacher,  and  was  found  to  be  1.9665  ±  0.0045  fm.   Rela- 
tivistic  corrections  proposed  by  Casper  and  Gross  [6,  18] 
were  compared  with  the  relativistic  correction  of  Friar  [8] 
and  no  significant  differences  were  found  in  the  region 
investigated. 
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I.   INTRODUCTION 

The  deuteron  being  composed  of  one  proton  and  one 

neutron  is  the  only  stable  two  body  nucleus  and,  as  such,  is 

the  simplest  of  the  nuclei  except  for  the  hydrogen  nucleus. 

The  proton  exists  as  a  stable  particle  in  the  form  of  a 

hydrogen  nucleus,  but  the  isolated  neutron  is  unstable.   In 

order  to  obtain  information  concerning  the  neutron,  this 

work  will  compare  electron  scattering  from  the  deuteron  and 

proton  and  extract  the  electric  form  factors  of  the  neutron, 

G„  ,  which  gives  information  on  the  charge  structure  of  the 

n 
neutron. 

2 

This  experiment  measured  G    at  various  values  of  q" 

n 
(the  square  of  the  Lorentz  invariant  four-momentum  trans- 
ferred from  the  scattered  electron  to  the  target  nucleus) 

-2  -2 

ranging  from  0 . 4  f m    to  0 . 5  f m   ,  in  order  to  resolve 

differences  among  several  earlier  experiments  [1,2,3,4]. 

The  results  of  this  work  augmented  with  those  of  Topping  [5] 

provide  values  of  the  neutron  charge  form  factor  for  the 

2-2-2 

momentum  range  q  =  0.0  5  fm    to  0.50  fm   .   By  comparing 

these  G„  ,  which  are  obtained  with  the  use  of  several  dif- 

n 
ferent  wave  functions,  to  the  G„   predicted  by  the  very 

n 
accurate  thermal  neutron  experiment  of  Krohn  and  Ringo  [4], 

a  model  that  best  represents  the  structure  of  the  deuteron 

may  be  found.   Three  different  models  containing  D-state 

contributions  of  4.57,  5.20  and  7.55%  were  tested  and 


although  the  results  are  not  conclusive,  the  data  favors 
the  model  with  7.55%  D-state.   Two  different  relativistic 
corrections  were  investigated,  the  first  developed  by 
Gross  [6]  and  modified  by  Schumacher  [7]  and  the  second 
developed  by  Friar  [8].   Results  using  both  methods  are 
reported  here.   For  the  momentum  transfers  of  this  experi- 
ment, it  is  shown  that  the  differences  between  the  two 

methods  are  negligible.   From  the  slopes  of  G^   for  the 

n 
various  deuteron  models  used,  a  mean  structure  radius  of 

the  deuteron  is  extracted  employing  a  method  developed  by 

Schumacher  [7]  . 

The  organization  of  this  thesis  is  as  follows:   the 
concepts  of  electron  scattering  and  the  definitions  of 
terms  needed  to  describe  the  scattering  from  protons  and 
deuterons  are  contained  in  Chapter  II. 

Chapter  III  discusses  briefly  the  experimental  equipment 
and  the  procedures  followed  in  order  to  obtain  -the  data. 

The  contents  of  Chapter  IV  discusses  the  methods  employed 
for  data  reduction.   Included  are  discussions  of  corrections 
and  normalizations  made  to  the  data,  radiative  and  ioniza- 
tion corrections  to  the  scattering  spectra,  and  a  detailed 
analysis  of  the  sources  of  errors. 

Chapter  V  presents  the  results  using  several  different 

wave  functions  and  both  methods  for  incorporating  the 

o 
relativistic  correction.   Plots  for  G    versus  q   are 

En 

included  for  all  cases.   Also  included  is  a  Rosenbluth  Plot 

as  an  experimental  verification  of  the  conventional  scaling 

law  assumptions  of  Chapter  II. 


Chapter  VI  presents  the  conclusions  of  this  work  and 
recommendations  by  the  authors  for  future  experiments.   Five 
appendices  amplify  details  of  several  topics  discussed  in 
this  thesis.   The  last  appendix  is  a  compilation  of  the 
FORTRAN  codes  used  in  the  data  analysis. 
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II.   THEORY 

A.   NEUTRON  CHARGE  FORM  FACTOR 

The  cross  section  for  an  extremely  relativistic  electron 
scattering  from  a  spinless  point  nucleus  with  no  magnetic 
moment  has  been  calculated  by  Mott  [9], 


Ve)  = 


r6a\ 


Mott 


Ze 

2E. 

l 


2  0 
cos   2   i 


.  4  e  s 

sin   j 


where 


r    ,    2Ei    .26 
%   =  1  +  — T  sin   -2- 

M  C 
n 


(II-D 


is  the  recoil  factor. 

The  experimental  cross  section  is  related  to  the  Mott 
cross  section  by 


d(A 


exp 


dcA 

d^Mott 


G2(q2)  , 


where  G (q  )  is  the  total  form  factor  of  the  nucleus,  which 
in  the  case  of  charge  scatterings  is  a  function  only  of  the 
square  of  the  four-momentum  transfer.   The  form  factor  takes 
into  account  the  finite  extent  of  the  nucleus. 

A  theoretical  expression  for  the  experimental  cross 
section  for  a  proton  has  been  developed  by  Rosenbluth  [10] 

and  is  given  in  terms  of  the  Mott  cross  section  and  total 

2 
form  factor  G   as  developed  by  Hand  et  al_.  [11] 
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§  =  aP,0,  G^,  (II-2, 

where  0^(8)  is  the  proton  Mott  cross  section  and 

2        2 

GE   +  T  GM 

G2  =  — E E  +  2t  tan2  |  G2  ", 

P       1  +  T  2   Mp 

where  G„  and  G,,   are  the  electric  and  magnetic  form  factors 
EM 
P       P  q2 

of  the  proton  respectively  and  x  =  -a— -5-  . 

4M 
P 

In  principle.  Equations  II-l  and  II-2  may  be  used  to 

measure  G    and  G   by  doing  experiments  at  constant  q  with 

P       P 
various  values  of  incident  energy  and  scattering  angle. 

However,  for  low  values  of  q,  the  magnetic  term  is  small  (at 
most  15%  of  the  total)  so  that  it  is  more  convenient  to  cal- 
culate the  magnetic  scattering  using  the  scaling  law  developed 
by  Lehmann  et  al.  [12]  and  others  [11],  which  was  shown  to 
be  valid  at  low  momentum  transfers  by  Brooks  and  Sessler  [13]. 
In  terms  of  the  proton  magnetic  moment,  u  /  the"  scaling  law 
is  given  by 

•   GM   =  ^p  GE   '  (II"3) 

P     r    P 

Utilizing  Equation  II-3,  Equation  II-2  can  be  simplified 
to 

4 

Gp  =  y£    [1  +  tu2  (1  +  2(1+t)  tan2  f}  ]  .       (II-4) 

2 
Since  the  tu   term  in  the  brackets  is  due  only  to  magnetic 

2  2  fi 

factors,  defining  C   =  tu  {1+2(1+T)tan   y}  ,  equation  II-4 

can  be  written 
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2    4 
gp  =  i+S  i**v  - 

Rewriting  this  and  defining  a  proton  magnetic  correction 
term,  Yp   =  . , -   ,  an  expression  for  G^   is  found, 


GE   -  Gp  Kl    '  (II"5) 

2 
From  Equation  II-2,  it  can  be  seen  that  G   can  be 

obtained  from  the  ratio  of  the  experimental  cross  section 

2 

to  the  Mott  cross  section  and,  G^   can  be  written  in  terms 

P 

of  measurable  quantities. 

doP 

GE   =  -f—   KM  (II"6) 

An  analogous  method  is  used  to  develop  an  expression  for 

2 

G^  ,  the  electric  form  factor  of  the  deuteron.   Jankus  has 

Ed 
developed  an  expression  for  the  cross  section  for  electron 

scattering  from  the  deuteron  [14].   Written  in  a  form  that 

includes  form  factors  to  account  for  the  finite  extent  of 

the  nucleons,  Jankus'  results  [15]  are 

d  - 2 


1+2  (1+n) tan      j 


da  „d,ft,  f_2    ,    2.       8      2_2     ,    2S    ,  2    _2     ,2. 

ST=    aM(6)lGE/q    )  +  9    n    GQ^(q    )  +3    V(q    } 
a  ad 

(II-7) 

2 

where  n  =  -3^-  and  G_  (q  )  ,  G_  (q  )  ,  GM  (q  )  are  respectively 

4Mj      Ed       Qd       Md 
a 

the  charge,  quadrupole  and  magnetic  moment  form  factors  of 
the  deuteron.   They  are  given  by: 

GE„  =  (GE   +  GE  '  CE  •    GQ„  "  <GE   +  GE  '  CQ  ' 
d      p     n  d      p     n 
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GM   =  2  ST  K   +  GE  »  CL  +  (GM   +  GM  '  CS  • 

d      p     p     n  p     n     ' 


Here,  C„  is  the  charge  structure  factor,  C_  is  the  quadru- 
pole  structure  factor,  CT  accounts  for  the  magnetic  contri- 
bution  to  scattering  arising  from  charge  convection  in  the 
deuteron  and  Cq  explains  the  contribution  to  scattering  by 
the  magnetic  moments  of  the  proton  and  neutron.   These  four 
factors  C^,  C^,  CT ,  and  C„  are  all  model  dependent  in  that 

h  \2  ±J  b 

they  arise  from  integrals  involving  S  and  D  state  wave 

functions  of  the  deuteron  [5],   Values  for  C„,  C_,  CT ,  and 

£•        U    1j 

Cq  for  momentum  transfers  studied  were  supplied  by  Lomon  [16] 
for  the  Lomon-Feshbach  Models  LF1,  LF5  and  LF15  [17]  which 
contain  percent  D-state  contributions  of  4.57,  5.20,  and 
7.55  respectively. 

If  u(r)  and  w(r)  are  the  normalized  S  and  D-state  wave 
functions,  then  Cp  is  given  by 


oo 

CE  =  /   [u2(r)  +  w2(r)]  jQ  (3p)  dr  , 


where  j   is  a  spherical  Bessel  function,  corresponding  expres- 
sions for  C  ",  C  ,  and    C     may  be  found  in  Ref.  [39].   In 
y   ij        o 

the  static  limit  (q=0),  the  above  functions  are  given  by 
GEd(0)  =  1  ,  GQd(0)  =  M2Q, 

GMd(0)=M^  CV%  "I  PD  VPn-T>]' 

cE(0)  =  1,  CQ(0)  =  m|q  , 

Cs(0)  =  1  -  3/2  PD  ,     CL(0)  =  3/2  PD  , 
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where  Q  is  the  quadrupole  moment  of  the  deuteron,  y   and  y 
are  the  proton  and  neutron  magnetic  moments  respectively, 

00 

2 
and  P   =  /  w  (r)dr  is  the  D-state  probability  of  the 

o 

deuteron. 

Relativistic  corrections  to  the  deuteron  scattering 
problem  have  been  calculated  by  Gross  [6] ,  and  Casper  and 
Gross  [18]  and  Friar  [8].   These  corrections  account  for 
relativistic  modifications  of  the  deuteron  wave  functions 
and  the  nucleon  current  incorporating  not  only  Lorentz 
contraction  in  two  dimensions,  but  also  the  fact  that  the 
equal-time  wave  function  in  the  rest  frame  is  not  the  same 
as  the  equal  time  wave  function  in  the  moving  frame. 

Schumacher  and  Bethe  [7]  and  Elias  et  al.  [15]  have  shown 

that  the  relativistic  corrections  of  Casper  and  Gross  [6,18] 

2 
to  the  deuteron  form  factors,  to  order  ■ °  ■  ,  enter  as  a 

2  4MP2 

factor  (1-  -t?)       on  the  right  hand  side  of    Equation  II-7. 

Friar  has  shown  that  the  relativistic  correction  enters 

in  two  ways;  first,  a  factor  (1+t)    should  multiply  the 

right  hand  side  of  Equation  II-7;  second,  C^  must  be  eval- 

2'     2 
uated  at  an  effective  q    =  q  /(1+n).   A  comparison  of  the 

results  using  both  formulations  may  be  found  in  Section  V. 

In  Appendix  D,  it  is  shown  that  by  expanding  Friar's 

2 
correction  in  a  binomial  series  for  low  q  ,  Schumacher's 

4 
correction  differs  by  less  than  4  parts  in  10   for  the 

2 

highest  q   investigated.   Consequently,  the  change  in  Cp  is 

2 

negligible  at  low  q  .   The  relativistic  corrections  of 

Casper  and  Gross  [18]  in  the  form  developed  by  Schumacher 
[7]  are  used  in  this  work. 

15 


When  Equation  II-7  is  used  to  analyze  the  deuteron 

experiments,  a  scaling  law  for  G    is  used  [3], 

Md 

Md 

V,  =  yd  M~  GE  , 
d       p    d 

-4 
The  quadrupole  term  is  dropped  since  it  is  less  than  5  x  10 

2 

in  comparison  to  the  total  at  the  highest  value  of  q 

reached  in  this  experiment  [38] .   Equation  II-7  may  be 
written  in  the  form 


an        -d   -  -  -a 


=   °M(6){GL   +   !  vl  4j[l+2U+n)    tan2   |]}(1    -   f)       . 


9  9   f~~  9     A  — 1 

Defining   C,   =   -_-  xy  ,  jl+2  (1+n )  tan      =■     and   noting    that 


do_d 
dft 


=   oM<6)    Gd    , 


where  G,  is  the  total  form  factor  of  the  deuteron,  then 

2 

Gd  =  (GE   +  GE  }  CE  (1  +  Cd)  U  "  l]        '  (H-8) 

P     n 

Additionally,  by  defining  KM   as  the  deuteron  magnetic  and 

r 
relativistic  correction  factor,  the  following  is  obtained: 

Kd 


Mr    1+Cd  \l-i/2j       ' 

From  Equation  II-8,  the  deuteron  electric  form  factor,  G„  , 

d  d 

is  obtained  in  terms  of  K„  ; 

M 
r 

GL    -    Gd-KM      ■    -fb  KM      =     <GE    +    GE    )2    CE     •  (II-9) 

d  r     .  aM(6)         r  p  n 
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The  ratio  of  the  proton  and  deuteron  electric  form  fac- 


tors, is  expressed  in  terms  of  G    and  G    as 

n       p 

6-  (ql)  (GL,  (q2)'+  Gj?     (qh  ) 


GE  V 
P   F 


JE  ^d* 


JE  %^d' 
n 


GE  <*p> 


(11-10) 


Combining  Equations  II-6,  and  II-9,  with  Equation  11-10, 
the  ratio  in  Equation  11-10  can  be  written  in  terms  of 
measured  quantities  only. 


tGE  (<Jd>  +  GE  "3d"  CE 
P 2 

GE  <<£ 

P     ^ 


da 


d  ^Hr- 


K 


an 

M    d  777 

r  GM<6) 


da 


P  -i 


K 


p  dQ 
Mc£(6) 
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It  is  important  to  note  that  q,  and  q  ,  which  are  the 
four  momenta  transferred  to  the  deuteron  and  proton  respec- 
tively have  different  values.   The  incident  energy  was 
adjusted  to  make  the  final  electron  energies  from  both  the 
deuteron  and  proton  scatterings  equal.   This  insures  that 
the  transmission  efficiency  of  the  spectrometer  and  the 
detector  efficiencies  were  constant  for  both  gases,  but  the 
four  momentum  transfers  were  slightly  different.   The 
deVries  b'  fit  [19]  for  the  proton  electric  form  factor  was 
used  to  adjust  for  this  slight  difference  in  momentum 
transfers . 
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The  ratio  in  Equation  11-11  may  be  written 
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Solving  this  for  G   (q,)  yields: 
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(11-13) 


R    is  the  experimental  ratio  of  the  deuteron  cross 
exp  c 

section   to  the  proton  cross  section,  and  R      is  the  ratio 

c  mott 

of  a„(8)  to  a£(6).   Note  that  af.(Q)    and  k£  were  calculated 
M         M  M  M 

at  the  incident  energy  of  the  electron  on  the  hydrogen,  which 
was  slightly  higher  than  the  energy  of  the  electrons  on  the 
deuterium. 
Writing 
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Equation  11-13  can  be  rewritten  in  the  more  convenient  form, 
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In  this  form  G^  (qj  includes  the  effects  of  the  Gross  rela- 
te   u 
n 

tivistic  corrections  to  the  deuteron  form  factor. 


B.   NEUTRON  ROOT  MEAN  SQUARED  CHARGE  RADIUS 

In  the  Breit  or  "brickwall"  frame,  the  form  factors 
G£   and  G   can  be  separated.   The  Breit  frame  is  the 


n 


n 
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reference  frame  in  which  the  incident  and  final  electron 
energies  are  equal,  and  the  initial  and  final  momenta  of 

the  deuteron  are  equal  but  oppositely  directed  [20]. 

2 

In  this  frame  G„  (q  )  is  the  Fourier  transform  of  the 

£i 

n 
charge  distribution  [21,22]  and  can  be  written  as 


oo  ,  -y     ■ 

(q2)  =  /   p(r)  eiq* 


GE  (q")  =  J  d3r 

n       o 


The  exponential  term  may  be  expanded  in  a  power  series 
about  q   =  0  and  the  coefficients  of  q   are  the  n   moments 
of  the  charge  distribution.   For  symmetric  charge  distribu- 
tions, only  the  even  moments  are  non-zero.   G    is  given  by 

n 

<r4    >    4 
GE    (q2)     =    GE     (0)     -    |   <r2R>    q2    +    J^jJL   +    ...f        (n-15) 
n  n 

2 

where  <r   >  is  the  mean  square  charge  radius  of  the  charge 
Cn 

2 
distribution.   The  slope  of  G    as  a  function  of  q   for  low 

Jit 

n 
momentum  transfers  is  related  to  the  mean  square  neutron 

charge  radius  by 

d   G_   -  -  i  <r?>  . 


,  2   E       6    CH 
dq    n 
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III.   EXPERIMENTAL  EQUIPMENT  AND  PROCEDURE 

A.  NAVAL  POSTGRADUATE  SCHOOL  ELECTRON  LINEAR  ACCELERATOR 
The  100  MeV  NPS  LINAC  is  an  electron  accelerator  similar 

in  design  to  the  Stanford  University  Mark  III  linear  accel- 
erator described  by  Chodorow  et  al.  [23]  .   The  thirty  foot 
accelerator  is  powered  by  three  Klystrons,  each  of  which 
receives  up  to  twenty-two  megawatts  of  peak  power.   The 
beam  energy  may  be  varied  from  30  to  105  MeV  with  an  average 
maximum  current  of  about  lyA.    The  machine  is  pulsed  at  a 
rate  of  60  pulses  per  second,  and  each  pulse  is  about  one 
microsecond  in  duration. 

The  beam  energy  is  selected  by  the  beam  deflection 
magnets  and  energy  defining  slits.   The  field  of  these  mag- 
nets is  measured  by  a  nuclear  magnetic  resonance  (NMR)  probe 
With  the  aid  of  quadrupole  magnets,  the  beam  is-  focused  to 
a  beam  diameter  of  about  2  mm  at  the  target.   A  detailed 
description  of  the  NPS  LINAC  is  given  by  Barnett  and 
Cunneen  [24] . 

B.  SPECTROMETER  AND  COUNTING  SYSTEM 

The  spectrometer,  described  in  detail  by  Oberdier  [25] , 
is  a  movable  120°,  16  inch  radius,  double-focusing  magnetic 
spectrometer  able  to  detect  electrons  scattered  at  angles 
between  45°  and  150°  in  steps  of  15°.   Ten  NE102  plastic 
scintillation  counters  span  part  of  the  focal  plane  of  the 
spectrometer,  where  each  detector  serves  as  the  source  of  a 
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counting  channel.   The  width  of  the  counters  allows  a  momen- 
tum spread,  Ap/p,  of  about  0.3%  for  each  channel.   Two  large 
detectors,  spanning  the  same  field  as  the  ten  smaller 
detectors,  served  as  dual  coincidence  backing  counters. 
Thus  an  electron  passing  through  both  backing  detectors  and 
a  front  detector  will  register  as  one  count  on  the  Jorway 
Model  1836  Dual  Channel  Scaler  located  in  the  control  room. 
A  triple  coincidence  pulse  will  register  on  a  scaler  only 
if  it  occurs  during  the  time  that  the  counter  gating  circuit, 
gated  with  the  accelerator  gun  pulse,  is  open.   This  pre- 
vents background  counts  from  registering  between  electron- 
beam  pulses. 

Scalers  are  capable  of  recording  periodic  pulses  up  to 
150  MHz,  but  due  to  the  pulsed  nature  of  the  beam,  the 
count  rate  was  held  below  60  counts/sec  on  all  runs  to  avoid 
large  count  rate  losses.   The  count  rate  correction  was 
found  to  be  linear  up  to  100  counts/sec  and  is  applied  in 
the  codes  BERHFORM  and  TRADFORM.   It  was  decided  to  use  only 
counters  4  through  7  in  order  to  consistently  follow 
Topping's  [5]  procedure  since  the  data  points  accumulated  by 
Topping  were  used  in  the  present  work.   Sufficient  data  were 
taken  such  that  the  ratios  found  from  channels  4  through  7 
were  well  defined. 

The  gas  target  cell  used  was  a  two  inch  diameter  steel 
cylinder  with  1.0  mil  Havar  windows  positioned  to  allow  the 
incident  electrons  to  enter  and  the  scattered  electrons  to 
exit.   Both  the  H_  and  D„  gases  were  maintained  at  liquid 
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nitrogen  temperature  and  about  150  psig.   The  gas  target 
system  is  described  in  detail  by  Topping  [5] . 

C.   DATA  ACCUMULATION 

Before  entering  into  the  actual  procedure  for  obtaining 
data  several  terms  must  be  defined.   A  "point"  is  the  set 
of  data  obtained  at  one  spectrometer  setting  for  either 
the  proton  or  deuteron.   It  consists  of  the  number  of 
counts  registered  in  each  of  the  ten  counters,  the  integra- 
tion capacitance  and  voltage/  machine  and  spectrometer 
energies  and  the  time  required  to  take  the  data.   A  "peak" 
is  the  collection  of  data  associated  with  one  of  the  ten 
channels,  e.g.,  "peak  4"  refers  to  the  data  collected  by 
counter  number  4.   A  "run"  is  the  collection  of  all  the 

points  of  all  the  peaks  for  both  the  deuterium  and  hydrogen 

2 

at  one  value  of  q   and  scattering  angle. 

2 

Topping  [5]  reported  results  of  21  runs  from  q   =  0.05 

-2  -2 

fm   to  0.35  fm   .   The  results  of  his  runs  and  the  12  runs 

made  by  the  authors  will  be  combined  to  extend  the  data  set 

-2 

to   .05  fm   .   The  data  collected  by  the  authors  were  mainly 

obtained  under  the  same  conditions  and  same  methods  as  used 
by  Topping  [5] . 

All  data  for  each  point  were  recorded  automatically  (see 
Section  IV-A) .   A  run  consisted  of  between  50  and  75  points 
first  taken  with  hydrogen  in  the  target  and  then  repeated 
with  deuterium. 

The  spectrometer  was  first  "saturated"  at  a  field  value 
corresponding  to  an  electron  energy  of  100  MeV  before  data 
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were  taken  on  either  H?  or  D„  to  insure  that  the  same 
hysteresis  loop  was  followed  on  both  spectra.   The  spectro- 
meter was  then  set  such  that  the  elastic  scattering  spectrum 
was  well  above  channel  4.   The  energy  selected  by  the 
spectrometer  was  stepped  down  0.0  5  MeV  between  each  point. 
This  process  was  continued  until  the  energy  selected  by 
the  spectrometer  was  about  two  full  widths  at  half  maximum 
below  peak  7.   Both  the  H^  and  D~  spectra  were  obtained  in 
the  same  manner. 

Two  additional  bits  of  data  were  recorded  by  hand,  the 
Secondary  Emission  Monitor  (SEM)  efficiency  and  the  pressure 
in  the  gas  target.   The  SEM  efficiency  was  measured  when 
each  of  the  peaks  registered  a  maximum  number  of  counts; 
therefore,  a  total  of  eight  values  were  taken  for  each  run 
and  they  varied  less  than  0.1%. 

The  above  data  were  always  taken  under  conditions  that 

were  as  similar  as  possible  for  the  deuteron  and  the  proton. 

2         -2  -2 

The  authors  made  runs  from  q   =  0 . 4  f m  '  to  0 . 5  f m   .   The 

time  required  to  complete  a  run  in  this  experiment  was 

greater  than  for  the  experiment  of  Topping  because  of  the 

2 

decrease  of  the  cross  section  with  q  ,  and  the  fact  that 

back  angles  had  to  be  used  exclusively.   Consequently  H_  and 
D2  spectra  were  accumulated  on  two  successive  days.   This 
was  a  deviation  from  Topping's  procedure,  but  this  method 
gave  data  that  reproduced  within  statistics. 

The  beam  charge  integration  was  always  set  such  that  the 
hydrogen  peaks  contained  from  10,000  to  20,000  counts  in  a 
single  channel. 
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As  in  Ref .  [5] ,  this  work  eliminates  the  need  for 

j    d  / 

counter  efficiencies  by  calculating  a  ratio  of  -jt^— / 

J  dft  /  dft 

for  each  channel. 
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IV.   DATA  ANALYSIS 

A.   DATA  REDUCTION 

The  data  collected  on  a  Model  33  teletype  was  transcribed 
onto  computer  data  cards  compatible  for  use  on  the  IBM  360-C. 
The  data  cards  were  organized   such  that  there  were  two  cards 
for  every  point  obtained  at  a  given  spectrometer  setting. 
The  first  card  contained  the  counts  for  each  of  the  ten  coun- 
ters, and  the  second  contained  the  number  of  background 
counts,  the  integration  time,  incident  electron  energy,  the 
voltage  integration,  the  energy  setting  for  the  spectrometer, 
the  capacitor  setting,  and  the  integrated  beam  charge  as 
measured  by  the  SFM.   The  hydrogen  and  deuterium  data  were 
treated  separately  with  FORTRAN  programs  BERHFORM.  and 
TRADFORM  respectively  which  are  listed  in  Appendix  E.   These 
programs  are  rivised  versions  of  TOP#FORM  [5].   They  are  used 
primarily  to  calculate  the  cross  section  measured  by  channels 
4  through  7,  but  they  also  make  count  rate  corrections  and 
charge  monitor  normalizations.   Charge  monitor  normalization 
was  required  because  the  automatic  termination  of  the  inte- 
gration for  each  of  the  spectrometer  settings  did  not  always 
stop  the  counting  system  at  the  same  value  of  beam  charge. 

The  correction  was 

V   , 
N (normalized)  =  N  (actual)   rererence  (IV-1) 

actual 

where  V      ,  is  the  actual  voltage  at  the  end  of  integration 
ac  uua  j_ 

and  V   ,.        is  an  arbitrary  standard  voltage  used  to  nor- 
rererence  J  3 

malize  the  hydrogen  and  deuterium  counts. 
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These  programs  were  also  used  to  plot  the  normalized 
count  rate  corrected  spectra  on  the  CALCOMP  plotters  at  the 
NPS  Computer  Center.   For  the  first  plotting  run,  certain 
additional  data,  which  are  necessary  to  calculate  cross 
sections,  (such  as  the  radiation  corrections  and  densities) 
were  given  arbitrary  values.   The  count-rate  corrected  and 
normalized  counts  for  both  deuterium  and  hydrogen  peaks 
were  plotted.   Using  these  plots,  the  position  of  each  peak 
and  the  line  shape  of  the  electron  spectrum  was  visually 
determined.   A  smooth  curve  was  fit  through  the  points  and 
those  points  that  were  affected  by  energy  instabilities  were 
appropriately  adjusted.   Typical  plots  of  both  raw  data  and 
smoothed  curves  are  reproduced  in  Figures  IV- 1,  IV— 2,  IV- 3 
and  IV- 4. 

The  radiation  correction  program  TOPRADCR  15]  requires 
the  following  inputs:   1)  the  peak  energies;  2)  the  pressure 
in  the  gas  target  for  hydrogen  and  deuterium;  and  3)  AE  which 
is  the  integration  interval  from  the  peak  energy  to  the 
final  spectrometer  setting  at  which  data  was  taken.   These 
data  were  obtained  from  the  before  mentioned  plots.   TOPRADCR 
also  used  the  composition  and  thickness  of  the  gas  target 
windows  and  the  virial  coefficients  of  hydrogen  and  deuterium 
as  stored  constants.   TOPRADCR  calculated  the  Schwinger 
radiation,  bremsstrahlung  and  Landau  ionization  corrections 
along  with  the  density  of  the  gas  for  each  channel  (peak) . 
The  necessity  for  the  radiation  corrections  will  be  discussed 
later  in  Section  IV-C.   Average  background  was  determined 
from  the  spectrum  above  the  elastic  peak. 
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The  total  radiation  and  ionization  corrections,  the 
gas  density,  the  SEM  efficiency  and  the  background  counts 
for  each  channel  were  typed  on  computer  cards  and  were  used 
as  input  data  to  BERHFORM  and  TRADFORM  along  with  the 
appropriate  data  decks  for  hydrogen  and  deuterium.   These 
programs  then  recalculated  the  elastic  cross  section  and 
the  statistical  error  in  that  cross  section  for  each  channel 
considered.   The  calculation  of  the  cross  section  is 
discussed  in  Section  IV-B,  and  the  method  for  calculating 
statistical  error  is  derived  in  Section  IV-D. 

Through  the  use  of  BERHFORM  and  TRADFORM,  four  cross 

sections  (one  for  each  channel)  were  obtained  for  each  of 

the  hydrogen  and  deuterium  peaks.   These  cross  sections 

were  then  put  in  the  form  of  a  ratio  (deuterium  cross 

section  divided  by  the  hydrogen  cross  section), and  a 

weighted  mean  ratio,  R  ,  was  calculated.   The  electric  form 
3  m 

factor  of  the  neutron,  G^  ,  for  each  run  was  determined 

n 

using  R  .   The  values  of  the  proton  electric  form  factors 
3   m  c 

and  the  proton  and  deuteron  Mott  cross  sections  and  magnetic 
corrections  were  obtained  from  programs  TOPROTON  and 
TOPDEUTR   [5] . 

B.   CALCULATION  OF  THE  ELASTIC  CROSS  SECTION 

The  experimental  cross  section  for  each  channel  is 

N    K   K   K.    -      _, 

°exp  "  N°  H-  AO  '  cm  ster   '  (IV-2) 

^     in   t 


31 


where   N    is  the  number  of  scattered  electrons  counted 
s  c 

in  a  particular  channel, 

N.   is  the  number  of  incident  electrons, 
in 

2 

N    is  the  number  of  target  atoms  per  cm  , 

Aft   is  the  solid  angle  defined  by  the  spectrometer 
slits , 

K  ,  K,  ,  and  K.  are  radiation  and  ionization  cor- 
s    b       1 

rections  (see  Section  IV-C)  , 
N    is  given  by 

pt  N 

N.  =      °  ,  (IV-3) 

t      A 

2 
where    p    is  the  density  of  the  target  (g/cm  ) , 

t    is  the  target  thickness  (cm) , 

N    is  Avogadro's  number,  and 

A    is  the  atomic  weight  of  the  target  gas. 

N.   is  determined  by 
in  J 

h.  -  ft-  _£3L_  ,  (iv-4) 

where    Q    is  the'  total  charge  of  incident  electrons, 

-19 

e    is  the  electron  charge  (1.60210  x  10     coulombs), 

V    is  the  voltage  of  the  integrator  in  volts, 

C    is  the  capacitance  of  the  integrator  in  farads , 

e    is  the  SEM  efficiency  (about  5.8%). 

sem  -1 

In  this  experiment  a  ratio  of  cross  sections  is  measured, 
so  that  several  of  the  parameters  of  Equations  IV- 2  and  IV-3 
which  are  not  absolutely  determined,  i.e.,  solid  angle  and 
target  thickness,  cancel  with  a  subsequent  elimination  of 
systematic  errors. 
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N    is  the  total  number  of  scattered  electrons, 
sc 

Area  under  scattering  peak 

„          (in  counts  -  MeV)  ,T.T  ._  > 

N    =  > ,  (IV-5) 

SC  6E 

where  E  is  the  energy  width  of  an  individual  detector. 

In  order  to  obtain  N   ,  the  number  of  counts  in  a  par- 
se e 

ticular  channel  must  be  integrated  over  the  energy  range 
of  interest. 

n-1 

N 


sc 


=      Z   i  [N.+N.+1]  (E.-E.+1)  /0. 0034877  E±,   (IV-6) 
i=l  ' 


where  N.  is  the  number  of  counts  in  the  channel  for  the  i 

1 

4-  Vi 

spectrometer  setting,  and  E.  is  the  energy  of  the  i    spec- 
trometer setting. 

Combining  Equations  IV-3,  IV- 4  and  IV-6  in  Equation  IV-2, 
the  experimental  cross  section  may  be  written  as 


n~1  i 1 

.  Z,    2  [Ni+Ni+l](Ei    Ei+1}  0. 00034  877E.  (1.0+0. 0039  (n- 5)) 
1=1  i 


s  b  l 


^P  CV  \    fpt   N 

o 


where  the  expression 


ee    /  Afi 

sem/  1   A   /  (IV-7) 


Ein  =  Ei  (1'°  +  °-0039  (n~5))  (IV-8) 


expresses  the  electron  energy  seen  in  counter  n  in  terms  of 
the  spectrometer  dial  setting  E..  Note  that  counter  5  sees 
the  electron  energy  read  from  the  spectrometer  control  dial. 

Count  rate  corrections  were  applied  to  N.  by  using  a 
linear  correction  empirically  determined  from  a  thick  target 
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experiment.   The  empirical  correction  is 


__  .  ,.       l  (recorded) 

N. (count  rate  corrected)  = 


'i^""1,  *-auc  wwJ.4.^w««,  -  1>Q  _0.00097  N./t.  ' 

1   1 

where  N./t.  is  the  average  count  rate.   Background  was  deter- 
mined  by  averaging  the  counts  for  each  peak  above  the  start 
of  the  elastic  peak.   Charge  integration  normalization  was 
done  as  given  in  Equation  IV-1. 

C.  is  defined  as  the  counts  in  the  i    channel/  corrected 
for  background  (B) ,  count  rate  and  charge  integration 
normalization  . 

N.  -^ 

i  vact 

Ci  =  1  -  0.00097  N./t  "  B  '  (IV_9) 

l 

The  final  expression  for  the  cross  section  in  terms  of 
defined  quantities  is 

,  n-1  (E.-E.  .  )  K  K,  K. 

i  i     (c.+c.Lj     1+1   s  b  x 


2  ._,  v  i   i+l;  0.0034877  (E .  /E .  ) 
exp  /  CV   \  /pt  N  \  '   UV  1U; 

-3 
where  Aft  is  the  solid  angle  and  has  a  value  of  1.831  x  10 

steradians . 

C.   RADIATION  CORRECTION  CALCULATIONS 

Both  the  proton  and  the  deuteron  cross  section  were 
derived  with  the  assumption  of  one  virtual  photon  exchange. 
This  process  is  represented  in  Figure  IV-5.   For  the  deuteron 
the  assumption  of  the  impulse  approximation  is  also  made,  as 
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Feynman    Diagram   of 
Rosenbluth    Calculation 

Figure  IV-5 


^ 

v 

•/x/^C         = 

1 

Xed  p 

<f 

1 

Feynman  Diagram   of 
Impulse  Approximation 

Figure  IV- 6 

X 

^ 
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represented  in  Figure  IV-6.   However,  in  the  actual  scatter- 
ing process,  more  processes  than  those  illustrated  in  the 
two  figures  occur.   The  experimenter  must  correct  for 
several  competing  processes  in  order  to  be  able  to  compare 
experiment  and  theory.   These  competing  processes  are 

1.  The  emission  of  real  photons  and  the  reabsorption 
of  virtual  photons  in  the  field  of  the  scattering  nucleus. 
The  correction  to  the  cross  section  for  this  process  was 
first  derived  by  Schwinger  [26]  and  is  therefore  called  the 
Schwinger  Correction. 

2.  Bremsstrahlung  emission  in  the  target  before  and 
after  the  nuclear  scattering  of  interest. 

3.  Excitation  and  ionization  of  material  in  the  target. 
In  order  to  obtain  the  cross  section  for  a  nuclear 

scattering  event  one  would  have  to  continue  down  to  zero 
energy  so  as  to  count  all  the  electrons  which  have  some 
probability  of  losing  all  their  energy  in  the  three  processes 
previously  discussed.   In  practice  the  measurement  of  the 
cross  section  is  stopped  at  some  AE  below  the  elastic  peak, 
so  it  is  necessary  to  include  a  correction  to  account  for 
the  undetected  electrons  with  energy  less  than  Ef  -  AE. 
Where  Ef  is  the  energy  of  the  elastic  peak  and  AE  is  the 
energy  from  the  elastic  peak  to  the  energy  corresponding  to 
the  last  spectrometer  setting,  usually  two  to  four  full 
widths  at  half  maximum. 
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C.   RADIATION  CORRECTION  CALCULATIONS 
1.   Schwinger  Correction 

The  Schwinger  correction  was  generally  about  15%  in 
this  experiment  and  is  the  largest  of  the  three  corrections. 

The  Schwinger  correction  accounts  for  processes 
which  occur  in  addition  to  the  one  virtual  photon  exchange. 
Schwinger' s  original  work  was  later  improved  by  Tsai  [27] 
and  by  Meister  and  Yennie  [28]  with  the  inclusion  of  the 
effects  of  proton  recoil.   The  correction  used  in  this  work 
was  evaluated  numerically  using  the  expression  taken  from  a 
review  article  by  Maximon  [29]. 
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=  exp 


2a 
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where 


E.  and  Er  are  the  electron  energies  before  and 
after  scattering, 

a  is  the  fine  structure  constant, 

m  is  the  electron  mass, 

2 

q   is  the  four  momentum  transfer  squared,  and 

/      n       /  n    •  \ 
dt  is  the  Spence  Function. 

0      t 

2 .   Bremsstrahlung  Correction 

When  an  electron  passes  through  material,  it  will 
radiate  energy  through  collisions  with  atomic  electrons  and 
nuclei  other  than  the  target  nucleus  of  interest.   This 
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energy  loss  which  is  referred  to  as  bremsstrahlung  has  been 
studied  by  Bethe  and  Heitler  [30].   The  energy  loss  per 
distance  of  material  traversed  is  given  by 

dx        o 

where  X   is  the  radiation  length,  which  is  the  mean  distance 
o  3 

over  which  the  electron  has  its  energy  reduced  by  a  factor 
of  1/e. 

Following  Butcher  and  Messel  [33]  the  radiation 
length  for  a  multi-element  target  was  calculated  using 

±J      i     )(znJL±i 

X  En.  A.    t   vi 

where  X   is  given  by 

A. 

XI  = 


°    4a  N   r2  Z(Z+0.8)  p.  £n(183  Z  1'3) 
o   e  l 

where   N   is  Avogadro's  number, 

r   is  the  classical  electron  radius, 

e  ' 

p.  is  the  density  in  g/cm  of  the  i    element,  and 
A.  is  in  grams /mole  of  the  itJl  element. 

Hofstadter  obtained  the  bremsstrahlung  correction 
factor  by  integrating  the  probability  function  P (a )  given 
by  Bethe  and  Heitler. 

(t/Jtn?) 
P(o)  =  _a , 

T  (l+t/£n2) 
where  P  (a)  is  the  probability  that  an  electron  traversing 
t  radiation  lengths  will  have  an  energy  loss  less  than  a 
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times  its  initial  energy  (aside  from  the  ionization  loss  and 
the  scattering  from  the  nucleus  of  interest) .   The  brems- 
strahlung  correction  factor  is 


Kfa  =  exp 


T    -   f  E 
£n 


X  £n2     lAE 


where    T  is  the  thickness  of  the  target  (this  includes 

the  two  steel  1  mil  windows  as  well  as  the  gas) 
in  g/cin2  ,  and 

E  is  the  geometric  mean  electron  energy  in  the 
target. 

In  this  experiment  the  bremsstrahlung  correction  is 
about  2%. 

3.   Ionization  Correction 

Electrons  also  lose  energy  by  excitation  and  ioniza- 
tion in  the  gas  and  in  the  gas  target  windows.  Landau  [34] 
calculated  the  average  energy  that  a  mono-energetic  beam 
would  lose,  and  the  spread  of  the  energy  distribution  of  the 
beam  after  passing  through  a  thin  target.  From. this  distri- 
bution it  is  possible  to  calculate  the  fraction  of  electrons 
which  are  not  detected  because  of  energy  degradation  through 
ionization. 

The  fraction,  i|>(AE),  of  the  beam  one  does  not  detect, 
because  of  ionization  losses  is 

E0  ,  E. 

r  2         J  r    mc 

ip(AE)  =  J    N(E)  dE/J      n(E)  dE  , 

o  /  o 

where  E„  is  the  energy  corresponding  to  the  energy  of  the 
last  spectrometer  setting, 

E.  _  is  the  incident  electron  energy,  and 
inc  ZIJ 
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N(E)  is  the  number  of  electrons  detected  at  energy  E 

Landau  [34]  derived  an  approximation  for  ip(AE)  ,  as 
have  others  [35],  which  is  used  in  this  work: 

iMAE)  -  1/u  / 
where  the  value  of  go  as  provided  by  Borsch-Supan  [36]  is 

£nA  +a\ 


o  *  ATI  - 


A  +  l 


where   A  =  y-1  =  -0.4228..., 
Y  is  Euler's  constant, 
A  =  AE/£  -  0.05 

|  n.  z. 

K   =    0.154  r      ^      *■    T  MeV,  and 

i  x   x 

2 

T  =  target  thickness  in  g/cm  , 

n .  is  the  number  of  atoms  of  the  i    element  per 
molecule  of  the  target,  and  Zj_  and  A^  are 
the  atomic  number  and  atomic  weight  of  the 
ith  element. 

The  ionization  correction  is  thus  given  by 
Ki  =  [1  -  MAE)]"1  • 

For  both  the  bremsstrahlung  and  ionization  correc- 
tions, the  target  was  treated  as  a  homogeneous  molecular 
compound  of  steel  and  either  H_  or  D„ ,  whereas  in  actuality 
the  target  consisted  of  H~  or  D^  gas  (2  inches  thick)  at 
77°K  and  150  psig  and  two  1  mil  steel  windows.   The  molecu- 
lar formulas  used  to  approximate  the  target  were  either 
Fe  H26  or  Fe  D„,  depending  on  the  gas  used. 
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D.   ERROR  ANALYSIS 

1.   Systematic  Errors 

In  this  experiment  the  proton  and  deuteron  cross 
sections  were  taken  with  the  same  geometry,  current,  and 
final  electron  energies.   At  most  a  time  delay  of  twelve 
hours  between  the  end  of  the  hydrogen  and  the  start  of  the 
deuterium  run  enabled  the  proton  and  deuteron  data  to  be 
taken  within  a  48  hour  time  span.   Consequently,  in  the 
ratio  of  the  deuteron  to  proton  cross  sections,  many  of  the 
systematic  errors  cancelled;  e.g.  solid  angle,  scattering 
angle,  SEM  efficiency,  target  thickness,  and  absolute 
spectrometer  efficiency. 

A  systematic  error  which  would  not  cancel  in  this  ratio 
experiment  is  that  which  arises  from  impurities  in  the 
hydrogen  and  deuterium  gases.   The  gases  used  in  this  experi- 
ment were  respectively  99.5%  deuterium  and  99.995%  hydrogen, 
as  determined  by  their  supplier  [37].   However,  the  equip- 
ment used  to  establish  these  purities  was  able  to  detect 
only  0.3%  or  more  contamination  of  D„  in  H_  or  H_  in  D~  . 
The  data  taken  in  this  experiment  used  the  same  bottle  of 
deuterium  and  two  different  bottles  of  hydrogen.   Consequently, 
the  hydrogen  data  may  have  varying  amounts  of  possible  con- 
taminants depending  on  which  bottle  was  used.   However,  if 
there  were  a  difference  in  the  contamination  of  H-  or  D„, 

this  difference  would  have  caused  a  non  zero  intercept  of 

2 
the  plot  of  G    versus  q  .   This  difference  was  not  evident 

n        2 

as  the  value  of  G^   (q   =  0)  is  zero  to  within  statistical 

n 
accuracy  for  all  deuteron  models  tested. 
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2 .   Random  Errors 

The  statistical  error  is  due  to  the  Poisson  distri- 
bution of  the  number  of  electrons  scattered  for  a  given 
number  of  incident  electrons.  This  error  was  calculated  by 
a  subroutine  in  the  BERHFORM  and  TRADFORM  programs  for  each 
individual  peak  and  included  the  statistical  errors  made  in 
determining  the  background  subtraction. 

If  X  is  the  cross  section,  the  following  may  be 
written 

M 

X  =  (Constant)  •   Z   (N   -  B)  , 

k=l    K 

where   N,  is  the  total  counts  (count  rate  and  voltage 

corrected)  in  the  ktn  spectrometer  setting, 

B  is  the  average  background  to  be  subtracted  and 
is  assumed  to  be  the  same  for  each  spec- 
trometer setting,  and 

M  is  the  number  of  spectrometer  settings. 

It  follows  that  if  N  =  N,  where  N  is  the  total  number  of 

k 

counts,  then: 

X  -  (Constant)  •  (N  -  M-B) 

The  variance  of  the  cross  section  is 

2  2 

2     8X     2     8X     2 

X     8N     N     9B    aB  ' 


a2  =  (Constant)2  (N  +  M2o?)     , 
a  a 

2 
where  a   is  the  variance  in  the  average  background,  and  is 

equal  to  B . .   Let  B   =  B.  (where  B.  is  equal  to  the  back- 
ground at  the  i^h  spectrometer  setting)  and  let  J  be  the 
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number  of  spectrometer  settings  above  the  peak  at  which 

J 

Z   B. 
.  _,   1 

background  was  taken.   Since  B  =  ,  then 


2     1   J    2      1   ?  n      1  „     B 

QB  =  72   *   V  =  3   *  Bi  =  72  BT  =  J  ' 
J   1=1    l    J   1=1       J 

which  in  turn  implies  that 

o  2  =  (Constant)2  (N  +  M2  ?r)  , 

x  J 

and  the  fractional  standard  deviation  of  the  cross  section  is 


QX 


N  -f  M2  B/J 
(N-MB) 2 


<h. 


Besides  the  statistical  (Poisson)  counting  errors,  the 
random  errors  in  this  experiment  originate  from  uncertainties 
in  determining  the  density  of  gas  in  the  target,  variations 
in  SEM  efficiency  during  a  run,  energy  shifts  of  the  inci- 
dent beam,  leakage  currents  in  the  beam  current  integrators, 
and  small  errors  in  the  radiation  corrections  resulting  from 
errors  in  determining  AE. 

Uncertainties  in  determining  the  density  of  the  gas  in 
the  target  were  caused  by  the  variation  of  pressure  in  the 
target,  that  is,  the  pressure  dropped  while  the  liquid 
nitrogen  was  being  pumped  and  then  rose  again  and  stayed 
steady  until  the  next  nitrogen  filling.   These  variations 
were  almost  identical  in  the  hydrogen  and  deuterium  runs 
and  were  about  0.3%  of  the  measured  pressure.   The  average 
pressure  for  each  peak  in  the  hydrogen  and  deuterium  runs 


43 


was  determined  separately  by  either  taking  the  pressure  at 
the  start  of  the  nitrogen  fill  or  at  the  end  of  the  fill 
depending  on  which  pressure  was  more  constant.   If  the 
pressure  at  the  start  (completion)  of  the  nitrogen  fill  was 
used  for  hydrogen,  then  the  pressure  at  the  start  (comple- 
tion) of  the  nitrogen  fill  was  used  for  deuterium.   The 
errors  in  the  measured  cross  sections  due  to  the  uncertainty 
in  determining  the  density  of  the  gas  in  the  target  varied 
from  run  to  run,  but  were  always  less  than  0.1%. 

The  efficiency  of  the  SEM  was  measured  four  times  during 
one  data  run  on  each  gas,  one  on  each  of  the  four  elastic 
peaks.   On  each  peak,  five  or  more  measurements  were  made. 
The  fractional  standard  deviation  of  a  single  measurement 
was  usually  less  than  0.2%,  giving  a  standard  deviation  of 
the  mean  of  the  SEM  readings  of  less  than  0.09%  for  the  four 
readings . 

Due  to  variations  in  line  voltage,  the  beam  energy  was 
not  constant,  and  as  a  result  the  scattered  energy  also 
changed  slightly  from  time  to  time  throughout  a  run.   Also 
there  was  a  slight  instability  in  the  spectrometer  power 
supply.   The  instabilities  had  the  effect  of  shifting  the 
counts  from  one  channel  to  another  on  the  counting  ladder. 
This  variation  in  the  data  was  most  noticeable  on  the  steep 
slopes  of  the  elastic  peaks  where  a  low  count  in  one  channel 
was  accompanied  by  a  high  count  in  the  adjacent  channel. 
It  was  found  that  by  closing  the  energy  defining  slits  more, 
the  variation  in  the  counts  due  to  shifts  in  energy  was 


44 


almost  negligible.   At  the  same  time,  it  was  noticed  that 
the  more  the  energy  defining  slits  were  closed  the  higher 
the  background  became.   Thus,  the  ratio  of  real  counts  to 
background  counts  became  less,  which  in  turn  increased  the 
errors  on  the  measured  cross  section.   The  position  for  the 
energy  defining  slits  was  adjusted  to  give  the  minimum 
combined  contribution  to  the  error  in  the  cross  section 
from  the  effects  of  energy  shifts  and  increased  background. 
Even  after  this  was  done,  in  some  cases  there  still  existed 
a  need  to  adjust  for  obvious  energy  shifts.   This  was 
accomplished  by  smoothing  the  data  points  to  the  general 
elastic  peak  line  shape.   It  is  estimated  that  in  the 
worst  case  these  energy  shifts  could  contribute  an  uncer- 
tainty in  the  average  measured  cross  section  of  as  much  as 
0.2%,  which  the  usual  error  in  the  measured  cross  section 
due  to  these  energy  shifts  was  about  0.1%. 

In  the  circuits  used  to  measure  and  integrate  the  beam 
current  there  existed  slight  leakage  currents.   If  these 
leakage  currents  were  not  constant  throughout  a  run,  then 
they  would  introduce  errors.   These  currents  were  measured 
and  found  to  be  constant  as  well  as  negligible  [5].   As  a 
result  the  only  error  attributable  to  these  leakage  currents 
existed  when  the  electron  beam  current  varied.   The  electron 
current  was  usually  held  constant  to  within  20%.   This  20% 
variation  would  cause  an  error  in  the  ratio  of  the  measured 
cross  sections  of  about  0.02%. 
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Another  source  of  uncertainty  is  the  inability  to 
determine  the  exact  position  of  the  peak  of  the  elastic 
spectrum  from  the  plotted  data.   This  uncertainty  is 
usually  about  0.01  MeV,  and  causes  an  error  in  the  Schwinger, 
bremsstrahlung,  and  ionization  corrections.   The  error  in 
the  measured  cross  section  due  to  this  AE  uncertainty  for 
all  three  radiation  corrections  was  of  the  order  of  0.05%. 
An  uncertainty  in  the  final  energy  at  which  one  terminates 
the  integration  of  the  elastic  spectrum,  contributes  approx- 
imately a  0.0  3%  error  to  the  measurement  of  the  cross 
section. 

A  list  of  all  the  errors  discussed  can  be  found  in 

Table  IV-1.   Excluding  statistical  error  (since  it  varies 

so  much  from  run  to  run) ,  the  average  error  in  the  cross 

section  ratio,  S^     _..    ,  is  about  0.24%.   The  total 

Avg.  Others 

error  of  the  weighted  mean  ratio  R  was  then  formed  from 

S?  .  ,  =  S2.  .  +  S2    _..  (IV-11) 

total     stat    Av.  Others 

One  significant  source  of  error  has  been  reduced  since 
Topping  [5]  performed  his  ratio  experiment.   It  was  found 
that  if  quadrupoles  3  and  4  were  used  instead  of  1,  2,  3  and 
4  that  the  background  was  effectively  50%  lower.   The  most 
plausible  explanation  for  this  reduction  of  background  is 
that  quadrupoles  1  and  2  defocused  the  electron  beam,  caus- 
ing the  beam  to  hit  the  beam  pipe  and  thereby  increased  the 
background.   The  data  in  this  experiment  was  collected  with 
only  quadrupoles  3  and  4  on,  and  as  a  result,  the  peak  to 
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TABLE  IV-1 
RANDOM  AND  SYSTEMATIC  ERRORS  IN  CROSS  SECTION  DATA 


Errors 


Average  Values  of  Errors 


Statistical 

Gas  Density 

SEM  Efficiency 

Energy  Shifts 

Leakage  Currents 

E  Measurement 
(Radiation  Corrections) 

E  Measurement 
(Integration) 

Average  Error 
(excluding  counting  statistics) 


0.3! 


0.1% 
0.09% 
0.  1% 
0.014% 
0.05% 

0.03% 


0.24,% 


background  ratio  of  counts  was  much  higher  than  in  the 
previous  work  [5] . 

The  uncertainty  of  G^   arises  from  the  errors  made  in 

■*      E 
n 

measuring  the  individual  cross  sections.   Let  R  =  o -,/a 

d'  p 

where  a   and  a ,  are  the  experimental  cross  sections  for  the 
P      d  e 

proton  and  deuteron  respectively.   Let  S   and  S,  be  the 

standard  deviation  of  the  proton  and  deuteron  cross  sections 

respectively.   G_,   is  determined  from 

n 
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^  d    p/   E 


-  1.0 


GE  (V   ' 
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(IV-12) 


and  G   /G   can  be  determined  from 
d   p 
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ge/ge 
d   p, 
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(IV-13) 


Solving  for  the  standard  deviation  in  R,  denoted  by  SR, 
due  to  the  errors  in  both  cross  sections,  one  obtains: 


(IV-14) 


The  standard  deviation  of  the  ratio  G   .  G   ,  denoted  by 

d    p 
Sn   /G  ,  is  obtained  from  Equation  IV-14, 
Gd   p 


S^  /G   = 

Gd   P 
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(IV-15) 


From  Equation  IV-12,  one  can  determine  the  standard 

deviation  in  G„  ,  denoted  by  S.,   ,  to  be: 

E  G,_, 
n  En 


E 


n 


CE      GE     R 


(IV-16) 


In  order  to  combine  the  results  obtained  from  four  different 

channels  in  one  run,  let  R.  denote  the  ratio  of  cross  sec- 

l 

tions  in  the  i    channel  and  let  SD   denote  the  standard 

th   1 
deviation  of  the  ratio  in  the  i    channel.   It  follows  that 

the  weighted  mean  ratio  R   can  be  calculated  from 

m 


R   = 
m 


7   f     2 

■\  K/sr. 

1=4  \ 
7 
Z        ll/S 


i=4 


Ri/ 


(IV-17) 
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where  the  standard  deviation  of  this  weighted  mean  ratio  is 
calculated  from 


R 


m 


L  E    l/s 


i=4 


R.  i 


(IV-18) 


R   and  S^   are  the  experimental  values  of  the  ratio  and 

m       R  e 

m 

of  the  standard  deviation  to  that  ratio. 
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V.   EXPERIMENTAL  RESULTS 

A.   NEUTRON  CHARGE  FORM  FACTORS 

The  21  runs  by  Topping  [5]  and  the  twelve  runs  by  the 
authors  are  reported  here.   Each  run  consisted  of  four 
measurements  of  R,  the  ratio  of  the  elastic  electron- 
deuteron  cross  section  to  the  elastic  electron-proton  cross 
section.   Of  these  132  measurements  of  R,  9  were  not 

included  due  to  reasons  given  in  Ref.  [5]. 

2 
Table  V-l  gives  the  q   transferred  to  the  deuteron, 

the  scattering  angle,  the  experimental  ratio  R  ,  and  the 

standard  deviation  of  R   determined  both  for  the  total 

m 

error  and  the  statistical  (Poisson)  error  for  a  given  run. 
Also  included  are  quantities  that  were  determined  by  cal- 
culations discussed  in  Section  II.   The  ratio  of  the 
deuteron  electric  form  factor  to  the  proton  electric  form 
factor,  derived  using  relativistic  corrections  by  both 
Friar  and  Gross,  is  also  listed  in  Table  V-l. 

Tables  V-2,3,4  and  5  present  the  neutron  charge  form 
factors  obtained  using  the  deuteron  models  LF1,5,  and  15 
with  the  relativistic  corrections  of  Gross  and  Friar.   The 
deuteron  structure  factors  used  are  presented  in  Appendix  D. 
The  best  estimates  of  the  error  in  each  of  these  33  values 

of  G    for  each  LF  model  and  relativistic  correction  is  the 

n 

error  obtained  using  S,_         ,  from  Equation  IV-11.   An  average 

r   total        ^  3 

2 
G^   and  the  associated  S,    ,  at  each  q   were  calculated 
E  total  M 

n 

using  Equations  V-l  and  V-2  respectively. 
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<GE  >. 

n  1 

7         2 
•     (S    ) 
(GE  )   =  i ^Ot_l  (v_1} 

n  m    £  x 

1  <Stot>i 


^to^m  = 


1  <Stot>iJ 


-% 


(V-2) 


2 

A  plot  of  the  average  Gp   (with  its  S    )  versus  q   is  pro- 

n 
vided  for  each  of  three  LF  models  and  both  relativistic 

corrections . 

All  values  of  G^   were  calculated  on  a  Wang  Model  700-C 

n 
calculator  using  programs  written  by  the  authors.   The  best 

2       2 

estimate  of  the  slope  of  G    versus  q  ,  d/dq   G   ,  was 

n  n 

obtained  from  a  linear  fit  weighted  by  S^_    for  all  values 

3      J      tot 

of  G„  .   For  comparison  purposes  the  linear  fit  was  also 

n 
forced  through  zero  for  each  plot. 

2 

The  average  value  of  G„   at  each  q  ,  using  Lomon- 

n 
Feshbach  Model  #1  is  given  in  Table  V-2. 

2 
The  average  value  of  G^   at  each  q  ,  usinq  Lomon- 

EN 

Feshbach  Model  #5  is  given  in  Table  V-3. 

2 

The  average  value  of  G„   at  each  q  ,    using  Lomon- 

h. 
n 

Feshbach  Model  #15  is  given  in  Table  V-4. 

The  value  of  G^   calculated  from  Lomon-Feshbach  Model 

n 
#15  and  its  SG    for  each  of  the  3  3  runs  is  given  in 

En 
Table  V-5. 
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B.  ROSENBLUTH  PLOT 

2  -2 

At  q   =  0 . 4  f m    deuteron  and  hydrogen  data  were  taken 

at  90°  and  120°  so  that  a  Rosenbluth  plot  could  be  con- 
structed.  The  theoretical  interpretation  of  the  Rosenbluth 
plot  is  discussed  in  Appendix  A.   Table  V-6  gives  the 
experimental  data  and  derived  quantities  used  in  the  Rosen- 
bluth plot  given  in  Figure  V-7.   The  use  of  the  scaling 
law  is  compatible  with  the  data. 

C.  CHARGE  RADIUS  OF  THE  DEUTERON 

Using  Method  II  given  by  Schumacher  and  Bethe  [7],  the 
deuteron  structure  rms  radius  was  calculated.   A  more 
detailed  explanation  is  given  in  Appendix  B.   The  structure 
radius  of  the  deuteron  for  LF1,  LF5,  and  LF15  was  provided 
by  Schumacher  and  Bethe  [7]  and  was  used  in  a  plot  of  the 


experimental 


T2   GE 
Ldq     nJ 


q2  =  0 


versus  r_,   (deuteron  structure 
Ed 


radius) .   This  plot  is  Figure  V-8  and  as  predicted  by 
Schumacher  such  a  plot  gives  a  straight  line  for  all  models 

with  the  correct  binding  energy.   The  intersection  of  this 

d 


line  with  the  value  of 


2  GE 

Ldq     nJ 


9    determined  by  the  very 

q^  =  0 

2 
accurate  neutron-electron  scattering  work  yields  r   .   The 

^d 

value  for  the  deuteron  structure  rms  radius  was  determined 

from  Figure  V-8  to  be  r^  =  1.9665  ±  0.0045. 

E 
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TABLE  V-l 

EXPERIMENTAL  RATIOS  AND  RESULTS  OF 
MODEL  INDEPENDENT  CALCULATIONS 


(fm  2) 


Exp.  Ratio 
d 

R  =  ^£ 
P 
exp 


Std.  Dev. 

of  R 
Total  and 
(stat. ) 


mott 


mott 


v# 


GE  «# 
P 


0.05 


0.10 


60° 


75' 


75 


75' 


75° 


0.9654 


0.9724 


0.9390 


0.9331 


0.9262 


0.0046 
(0.0040) 

0.0038 
(0.0030) 

0.0042 
(0.0035) 

0.0034 
(0.0024) 

0.0046 
(0.0038) 


0.30269 
0.31356 

0.16837 
0.17578 

0.81525 
0.86757 


0.99400 


0.99396 


0.98789 


0.81525J  0.98789 
0.867571 


0.81525 
0.86757 


0.98789 


0.99404 


0.99404 


0.98815 


0.98815 


0.98815 


2 
^d 

(fm-2) 


K   (Gross) 
r 


K 


K   (Friar) 
r 


M 


E 


(Gross) 

Total 

and 
(stat) 
Error 


+ 


(Friar) 


Total 
and 
GE   -(stat) 
p    Error 


+ 


0.05 


0.10 


60 


75° 


75° 


75' 


75' 


1.00001 
1.00001 

0.99997 
0.99996 

0.99993 
0.99992 

0.99993 
0.99992 

0.99993 
0.99992 


0.99338 
0.99114 
0.98230 
0.98230 
0.98230 


0.9685±0.0023 
(0.0020) 

0.9693+0.0019 
(0.0015) 

0.9475±0.0021 
(0.0018) 

0.9445±0.0017 
(0.0012) 

0.9410±0.0023 
(0.0019) 


0.9685±0.0023 
(0.0020) 

0.9693±0.0019 
(0.0015) 

0.9475±0.0021 
(0.0018) 

0.9445±0.0017 
(0.0012) 

0.9410±0.0023 
(0.0019) 
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TABLE  V-l   (Continued) 


(fm  2) 


Exp.  Ratio 

d 
o 

r  =  -£2B 

exp 


Std.  Dev. 

of  R 
Total  and 
(Stat. ) 


mott 


mott 


GE  «Jp> 

P 


GE  (*d> 

P 


0.10 


0.20 


75' 


60 


60 


75° 


90 


90' 


0.9221 


0.8204 


0.8402 


0.8323 


0. 8365 


0.0040 
(0.0031) 

0.0051 
( o  .  n  o  4  5 ) 

0.0076 
(0.0072 

0.0049 
(0.0042) 

0.0049 
(0.0042) 

0.0052 
(0.0046) 


0.81525 
0.86757 

0.70815 
0.76120 

0.70815 
0.76120 

0.39041 
0.42650 

0.22434 
0.24871 

0.22434 
0.24871 


0.9878Q 
0.97595 

0.97595 
0.9758? 
0. 97569 
0.07569 


0.98815 


0.97653 


0.Q7653 


0.97653 


0.97653 


0.97653 


K^    (Gross) 

r 


(fm"2) 


K 


KT.    (Friar) 
Mr 


M 


(Gross) 

Total 

jd   +    and 
"  (Stat) 
Jp      Error 


(Friar) 


Total 

and 
(Stat) 
Error 


0.10 


0.20 


75° 


60° 


60' 


75' 


90 


90' 


0.99993 
0.99992 

1.00001 
1.00001 

1.00001 
1.00001 

0.99985 
0.99985 

0.99896 
0.99896 

0.99896 
0.99896 


0.98230 


0.97359 


0.97359 


0.96493 


0.9512Q 


0.95129 


0.9389±0.0020 
(0.0016) 

0.884810.0028 
(0.0024) 

0.888710.0041 
(0.0039) 

0.892110.0026 
(0.0022) 

0.887110.0026 
(0.0022) 

0 . 889310.0028 
(0.0024) 


0.938910.0020 
(0.0016) 

0.884810.0028 
(0.0024) 

0.888710.0041 
(0.0039) 

0.892110.0026 
(0.0022) 

0.887110.0026 
(0.0022) 

0.  88°310.0028 
(0.0024) 
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TABLE  V-l   (Continued) 


2 

^d 

(fm-2) 


Exp.    Ratio 

d 
a 

R  =  -^E 

ap 
exp 


Std.    Dev. 

of    R 
Total    and 
(Stat. ) 


HiOtt 


mott 


GE     <1p> 

P 


GE     «# 
P 


0.20 


0.25 


0.30 


120° 


120' 


120° 


75° 


75° 


90 


0.8068 


0.8280 


0.8057 


0.7987 


0.6973 


0.7718 


0.0118 
(0.0106) 

0.0102 
(0.0099) 

0.0082 
(0.0079) 

0.0056 
(0.0051) 

0.0049 
(0.0043) 

0.0057 
(0.0052) 


0.72058 
0. 81829 

0.71893 
0. 81637 

0 .71893 
0. 81637 

0. 30690 
0. 33890 

0.27183 
0.28073 

0.14386 
0. 16338 


0.97555 


0.97550 


0.97550 


0.96982 


0.96625 


0.96365 


0.97653 


0.97653 


0.97653 


0.97082 


0.96516 


0.96516 


2 

^d 

(fm'2) 


K^    (Gross) 
r 


(Gross) 


(Friar) 


K 


1C    (Friar) 


M 


Total 

and 
(Stat) 
Error 


E 


Total 

and 
(Stat) 
Error 


0.20 


0.25 


0.30 


120c 


120c 


12  0C 


75' 


75° 


90 


0.99466 
0.99466 

0.99465 
0.99464 

0.99465 
0.99464 

0.99982 
0.99981 

0.99978 
0.99977 

0.99845 
0.99844 


0.89020 


0.88999 


0. 88999 


0.95641 


0.95135 


0.92817 


0.8901±n.0065 
(0.0058) 

0.9018±0.0055 
(0.0054) 

0.8896±0.0O45 
(0.0044) 

0.8687+0.0030 
(0.0028) 

n.8426±0.0O3O 
(0.0026) 

0.8538±0.0031 
(0.0029) 


0.8901+0.0065 
(0.0058) 

0.9018+0.0055 
(0.0054) 

0.8895+0.0045 
(0.0044) 

O.8687±0.O03O 
(0.0028) 

0. R426±0.n030 
(0.0026) 

0.853  8+-  0.0031 
(0.002°) 
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TABLE  V-l   (Continued) 


(fm  2) 


Exp.  Ratio 

d 
a 

r  =  -S2E 

exp 


Std.  Dev. 

of  R 
Total  and 
(Stat. ) 


mott 


mott 


GE  <1p> 
P 


P 


0.30 


0.35 


0.40 


90 


90 


90° 


90' 


90° 


90' 


0.7529 
0.7558 
0.7192 
0.6596 
0.6815 
0.6829 


0.0060 
(0.0055) 

0.0070 
(0.0066) 

0.0054 
(0.0048) 

0.0056 
(0.0051 

0.0056 
(0.0050 

0.0058 
(0.0053) 


0.14386 
0.16338 

0.14386 
0.16338 

0.12122 
0.13916 

0.13341 
0.13916 

1.04376 
1.210236 

1.04376 
1.21024 


0.96365 

0.96365 
0.95767 
0.96112 
0.9  51718 
0.95172 


0.96516 
0.96516 
0.95955 
0.95955 
0.95400 
0.95400 


2 
<*d 

(fm-2) 


K„  (Gross) 

M 

r 


(Gross) 


K 


KM  (Friar) 
Mr 


M 


Total 

and 
(Stat) 
Error 


(Friar) 


Total 

and 
(Stat) 
Error 


0.30 


0.35 


0.40 


90° 


90' 


90 


90' 


90 


90' 


0.99845 
0.99844 

0.99845 
0.99844 

0.99819 
0.99818 

0.99819 
0.99818 

0.997933 
0.997918 

0.997933 
0.997918 


0.92817 


0.Q2817 


0.91695 


0.92340 


0 .  905924 


0.905^24 


0.8433±0.0O34 
(0.0031) 

n. 844°±0.0039 
(0.0037) 

0.8242±0.0031 
(0.0027) 

0.828110.0035 
(0.0032) 

0.802710.0033 
(0.0030) 

0.803610.0034 
(0.0031) 


0.843310.0034 
(0.0031) 

0. 844910.0039 
(0.0037) 

0.824210.0031 
(0.0027) 

0.828110.0035 
(0.0032) 

0.802710.0033 
(0.0030) 

0.803610.0034 
(0.0031) 
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TABLE   V-l       (Continued) 


(fm~2) 


Exp.    Ratio 

d 

a 

R  =  -£2E 

aP 
exp 


Std.    Dev. 

of    R 
Total    and 
(Stat. ) 


mott 


mott 


GE     «Jp> 

P 


GE     «£ 
P 


0.40 


0.477 


120c 


12  0° 


120c 


90' 


90° 


90' 


0.6430 


0.6266 


0.6300 


0.6324 


0.6389 


0.6330 


0.0062 
(0.0057) 

0.0061 
(0.0056) 

0.0062 
(0.0057) 

0.0061 
(0.0056) 

0.0051 
(0.0045) 

0.0060 
(0.0055) 


0.32879 
0.39489 

0.32879 
0.39489 

0.32879 

0.3Q48Q 

0.85522 
1.00679 

0.85522 
1.00679 

0.85522 
|1.00679 


0.95117 
0.95117 
0.95117 
0.94262 
0.94262 
0.94022 


0.95400 


0.95400 


0.95400 


0.94550 


0.94550 


0.94500 


K^J    (Gross) 
r 


(Gross) 


(fnf2) 


K 


KM    (Friar) 


M 


Total 

and 
(Stat) 
Error 


(Friar) 


Total 

and 
(Stat) 
Error 


0.40 


0.477 


120° 

120° 

120° 

90° 

90° 

90° 


0.98938 
0.98936 

0.98938 
0.98936 

0.98938 
0.98936 

0.99754 
0.99752 

0.99754 
0.99752 

0.99754 


0.798^4 


0.79894 


0.79894 


0.88935 


0.88935 


0.88935 


0.8118±0.0039 
(0.0036) 

0.801410.0^39 
(0.0036) 

0.8036±0.0039 
(0.0036) 

0.7739±0.0037 


0.7778±0.0031 
(0.0027) 

0.7742±0.0037 
(0.0034) 


0.8118±0.0039 
(0.0036) 

0. 8014±0.0039 
(0.0036) 

0.8036±O.OO39 
(0.0036) 

0.7739±0.0037 
(0.0034) 

0.7778±0.0031 
(0.0027) 

0.7742±0.0O37 
(0.0034) 
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TABLE   V-l       (Continued) 


(fm~2) 


Exp.     Ratio 

d 

a 

R  =  -S2E 

a? 
exp 


Std.    Dev. 

of    R 
Total   and 
(Stat. ) 


mott 


mott 


GE     <"p> 
P 


P 


0.50 


105c 


105' 


10  5° 


105< 


0.6033 


0.6042 


0.6044 


0.6015 


0.0057 
(0.0052) 

0.0057 
(0.0052) 

0.0058 
(0.0053) 

0.0060 
(0.0055) 


0.46122 
0.55662 

0.46122 
0.55662 

0.46122 
0.55662 

0.46122 
0.55662 


0.93952 


0.93952 


0.93^52 


0.93952 


0.94307 


0.94307 


0.94307 


0.94307 


K^    (Gross) 


(fm~2) 


K 


K       (Friar) 

Mr 


M 


(Gross) 


E 


Total 

and 
(Stat) 
Error 


(Friar) 


Total 

and 
(Stat) 
Error 


0.50 


105° 


105° 


105° 


105' 


0.99368 
0.9  9366 

0.99368 
0.99366 

0.99368 
0.99366 

0.99368 
0.99366 


0.83642 


0.  83642 


0.83642 


0.83642 


0.7677+0.0036 
(0.0033) 

0.7683±0.0036 
(0.0033) 

0.7685±0.0037 
(0.0034) 

0.7666+0.0938 
(0.0035) 


0.7677±0.0036 
(n.0033) 

n.7683±0.0036 
(0.0033) 

0.76R5±0.0037 
(0.0034) 

0.7666±0.0038 
(0.0035) 
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TABLE  V-2 

AVERAGE  G    AT  EACH  q2  USING  LOMON-FESHBACH  MODEL  #1 
n 


2 

q 

Gross*  Relativistic 
Correction 

Friar's  Relativistic 
Correction 

(fnf2) 

GE 
n 

GE 
n 

GEn 

0.05 

0.0002 

0.0015 

0.0002 

0.0015 

0.10 

0.0037 

0.0011 

0.0037 

0.0011 

0.20 

0.0027 

0.0013 

0.0027 

0.0013 

0.25 

0.0070 

0.0034 

0.0069 

0.0034 

0.30 

0.0076 

0.0019 

0.0075 

0.0019 

0.35 

0.0093 

0.0027 

0.0092 

0.0027 

0.40 

0.0084 

0.0020 

0.0081 

0.0020 

0.477 

0.0098 

0.0025 

0.0095 

0.0025 

0.50 

0.0106 

0.0023 

0.0102 

0.0023 

A  linear  fit  to  the  data  in  Table  V-2  gives: 
G„ 


n 


n 


=  (0.0003  ±  0.0010)  +  (0.0208  ±  0.0040)q' 
for  Gross1  correction  and 

=  (0.0004  ±  0.0010)  +  (0.0200  ±  0.0040)q' 
for  Friar's  correction. 


A  linear  fit  of  the  same  data  forced  through  zero  gives: 

GE  (q2)  =  (0.0218  ±  0.0022)q2  for  Gross1  correction  and 
n 

G   (q  )  =  (0.0213  ±  0.0022)q2  for  Friar's  correction. 


E 


n 
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TABLE  V-3 


AVERAGE  G^   AT  EACH  q   USING  LOMON-FESHBACH  MODEL  #5 


n 


2 

q 

Gross ' 

Relativistic 

Friar  *  s 

Rela 

tivistic 

Corre 

ction 

Correct 

ion 

(fm-2) 

GE 

n 

Sge 

GE 

n 

S°E 

^n 

0.05 

0.0002 

0.0015 

0.0002 

0.0015 

0.10 

0.0037 

0.0011 

0.0037 

0.0011 

0.20 

0.0028 

0.0013 

0.0027 

0.0013 

0.25 

0.0070 

0.0034 

0.0069 

0.0034 

0.30 

0.0076 

0.0019 

0.0075 

0.0019 

0.35 

0.0094 

0.0027 

0.0092 

0.0027 

0.40 

0.0084 

0.0020 

0.0082 

0.0020 

0.477 

0.0099 

0.0025 

0.0096 

0.0025 

0.50 

0.0107 

0.0023 

0.0103 

0.0023 

A  linear  fit  to  the  data  in  Table  V-3  gives: 
'E 


G„  (q2)  =  (0.0003  ±  0.0010)  +  (0.0210  ±  0.0040)q2 


n 


for  Gross'  correction  and 


G   (q2)  =  (0.0004  ±  0.0010)  +  (0.0202  ±  0.0040)q2 


E 


n 


for  Friar's  correction. 


A  linear  fit  of  the  same  data  forced  through  zero  gives: 

G   (q  )  =  (0.0220  ±  0.0022)q2  for  Gross'  correction  and 
n 

G   (q  )  =  (0.0214  ±  0.0022)q2  for  Friar's  correction. 


E 


n 
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TABLE  V-4 

AVERAGE  G^   AT.  EACH  q2  USING  LOMON-FESHBACH  MODEL  #15 
n 


2 

q 

Gross*  Relativistic 
Correction 

Friar's  Relativistic 
Correction 

(fm"2) 

GE 
n 

^n 

GE 

n 

En 

0.05 

0.0000 

0.0015 

-0.0000 

0.0015 

0.10 

0.0033 

0.0011 

0.0033 

0.0011 

0.20 

0.0020 

0.0013 

0.0020 

0.0013 

0.25 

0.0061 

0.0034 

0.0060 

0.0034 

0.30 

0.0066 

0.0019 

0.0066 

0.0019 

0.35 

0.0082 

0.0027 

0.0080 

0.0027 

0.40 

0.0070 

0.0020 

0.0068 

0.0020 

0.477 

0.0083 

0.0025 

0.0080 

0.0025 

0.50 

0.0090 

0.0023 

0.0086 

0.0023 

A  linear  fit  to  the  data  in  Table  V-4  gives: 

G,,  (q2)  =  (0.0002  ±  0.0010)  +  (0.0177  ±  0.0040)q2 


E 


n 


for  Gross '  correction  and 


Gp,  (q2)  =  (0.0003  ±  0.0010)  +  (0.0170  ±  0.0040)q2 


n 


for  Friar's  correction. 


A  linear  fit  of  the  same  data  forced  through  zero  gives 
;E 


2  2 

G   (q  )  =  (0.0185  ±  0.0022)q   for  Gross'  correction  and 


n 


G   (q2)  =  (0.0180)  ±  0.0022)q2  for  Friar's  correction. 


n 
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TABLE  V-5 
FOR  EACH  OF  THE  30  RUNS 
CALCULATED  FROM  LOMON-FESHBACH  MODEL  #15 


G„   AND  S.  . 
E        tot 
n 


2 

q 

Gross ' 

Relativistic 

Friar ' s 

Relativistic 

Correction 

Correction 

(fm-2) 

GE 

S% 

GE 

Sr 

n 

En 

n 

^n 

0.05 

0.0003 

0.0019 

0.0003 

0.0019 

-0.0005 

0.0024 

-0.0005 

0.0024 

0.10 

0.0078 

0.0022 

0.0078 

0.0022 

0.0010 

0.0025 

0.0010 

0.0025 

0. 0047 

0.0018 

0.0047 

0.0018 

-0.0012 

0.0021 

-0.0012 

0.0021 

0.20 

-0.0027 

0.0030 

-0.0028 

0.0030 

0. 0015 

0.0045 

0.0014 

0.0045 

0.0052 

0.0029 

0.0052 

0.0029 

-0.0002 

0.0029 

-0.0003 

0.0029 

0.0021 

0.0030 

0.0021 

0.0030 

0.0030 

0.0072 

0.0030 

0.0072 

0.0159 

0.0061 

0.0158 

0.0061 

0.0025 

0.0050 

0.0024 

0.0050 

0.25 

0.0061 

0.0034 

0.0060 

0.0034 

0.30 

0.0153 

0.0036 

0.0151 

0.0036 

0.0032 

0.0039 

0.0030 

0.0039 

0.0050 

0.0045 

0.0049 

0.0045 

0.0024 

0.0034 

0.0023 

0.0034 

0.35 

0.0061 

0.0036 

0.0060 

0.0036 

0.0108 

0.0041 

0.0106 

0.0041 

0.40 

0.0050 

0.0039 

0.0048 

0.0039 

0.0060 

0.0041 

0.0058 

0.0041 

0.0159 

0.0047 

0.0157 

0.0047 

0.0034 

0.0046 

0.0032 

0.0046 

0.0060 

0.0047 

0.0058 

0.0047 

0.477 

0.0061 

0.0046 

0.0058 

0.0046 

0.0110 

0.0038 

0.0107 

0.0038 

0.0066 

0.0045 

0.0063 

0.0045 

0.50 

0.0089 

0.0045 

0.0086 

0.0045 

0.0096 

0.0045 

0.0092 

0.0045 

0.0098 

0.0046 

0.0095 

0.0046 

0.0074 

0.0047 

0.0071 

0.0047 
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Figure  V-8 
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VI.   CONCLUSIONS 

As  a  result  of  the  123  measurements  of  the  ratio  of  the 
elastic  electron-deuteron  cross  section  to  the  elastic 

electron-proton  cross  section,  nine  values  of  G„   at  dif- 

2  n 

ferent  values  of  q   were  determined  for  each  deuteron  wave 

2 

function  considered.   From  these  values  of  G^   and  q  ,  a 

n 


7    GF 

.dq     n 


[; 


„     was  determined. 
q^  =  0 


Relativistic  corrections  proposed  by  both  Friar  [8]  and 
Gross  [6]  were  investigated.   It  was  found  that,  for  a  given 


deuteron  wave  function,   ^   G 


2   E 
dq     n  I 


„    was  smaller  using  the 


q  =0 
Friar  correction.   Although  this  pattern  was  consistent,  it 

was  also  found  that  the  difference  in  the  two  slopes  was 

always  well  within  one  standard  deviation.   No  attempt  was 

made  to  select  one  form  of  the  relativistic  correction  as 

best  because  the  difference  was  so  small. 

All  wave  functions  considered  (LF1,  5,  15)  were  in 

agreement  with  the  very  accurate  thermal  neutron  slope  of 


(0.0193    .0004)  fm2.   The  values  of   -^r-  GE 

,2    n 
•dq 


4 


2  _  obtained 
q  =0 


using  Gross1  relativistic  corrections  were:   for  LF1  (0.0218  ± 
±  0.0022),  for  LF5  (0.0220  ±  0.0022)  and  for  LF15  (0.0185 
±  0.0022).   The  values  obtained  using  Friar's  corrections 
were:   for  LF1  (0.0213  ±  0.0022),  for  LF5  (0.0214  ±  0.0022), 
and  for  LF15  (0.0180  ±  0.0022).   LF15  was  found  to  be  in 
closest  agreement  with  the  thermal  neutron  slope. 
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A  value  for  the  structure  radius  of  the  deuteron  was 
also  obtained  from  this  data.  This  value  was  determined 
by  comparing  the  slopes  given  by  each  wave  function.   The 

structure  radius  was  determined  to  be  r_  =  (1.9665  ±  0.0045)fm, 

2  -2 

A  Rosenbluth  plot  was  made  at  q   =  0 . 4  f m   .   The  results 

of  this  plot  are  discussed  in  Appendix  A  and  show  agreement 
with  the  use  of  the  conventional  deuteron  scaling  law. 

If  this  experiment  were  to  be  repeated  or  extended  it  is 

2 

recommended  that  data  at  higher  q   be  taken  in  order  to 

4 
determine  more  accurately  the  coefficient  of  the  q   term  in 

the  expansion  of  G^  .   In  Appendix  C  a  coefficient  for  the 

hi 

4  n 

q   term  was  taken  from  work  done  at  DESY  [38]  and  it  was 

2  -2 

found  that  at  q   =  0  .  5  f  m    this  term  decreases  G^   by 

n 

about  10%.   The  effect  on  the  slope  of  the  G„   curve  is  also 

n 

shown  to  be  about  a  7%  increase. 


73 


APPENDIX  A 
THE  ROSENBLUTH  PLOT 

Based  on  the  one  photon  exchange,  the  cross  section  for 
electron-proton  scattering  was  derived  by  Rosenbluth  [10]. 
Using  the  same  assumption  (one  photon  exchange)  Jankus  [14] 
calculated  the  deuteron  cross  section.   Jankus1  result   is 
given  in  Equation  II-7  and  is  reproduced  below 


da 


da 


d 


dU 


dft 


fmott 


_2  ,  8   2  _2     2    r2 

GE/  9  n  G0.    +  3  n  GM, 
d  d  d  L 


1+2  (1+n)  tan 


2  0 


where  G-,  can  be  written  as 
d 


<- 


,-2,2,    8   2  „2     2    _2 

GE <*  '  +  9  1   V  +  3  n  GM„ 

d  ^d         d 


2  0*1 
1+2  (1+n)  tan   §■  I 

2  J 


Since  the  deuteron  cross  section  can  be  measured  at  different 
scattering  angles  while  at  the  same  four  momentum  transfer, 

one  can  see  that  a  plot  of  G,  versus  (1  +  2  (1+n)  tan   ■») 

2     2 

should  give  a  straight  line  with  -^  n  G    as  the  slope  and 

2     8   2   2 

G    +  g-  n   G    as  the  intercept.   By  plotting  the  deuteron 

d         yd 
total  form  factor  G,  versus  [1+2  (1+n) tan   y]  ,  the  adequacy 

of  the  conventional  scaling  law  can  be  tested.   Such  a  plot 

2  -2 

has  been  made  for  q   =  0 . 4  f m    and  it  can  be  found  in 

Section  V. 

One  cannot  use  Equation  II-7  directly  since  it  contains 

the  absolute  cross  section  of  the  deuteron.   However  the 

deuteron  cross  section  can  be  obtained  by  using  the  hydrogen 
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cross  section  along  with  the  deVries  b'  fit.   The  deVries  b' 
cit  [19]  is  used  because  it  is  the  best  fit  to  the  absolute 
proton  cross  section.   It  is  used  as  a  normalization  factor 
with  the  experimental  cross  section  obtained  from  the  data. 
That  is, 


doJ 
dfi 


N 


daJ 


'deVries 


dft 


fexp 


Therefore, 


N  = 


daJ 


daJ 


dft 


deVries 


dfi 


exp 


N  is  a  normalization  factor  for  a  given  experiment  at  a 
given  angle  and  energy.   Since  both  the  proton  cross  section 
and  the  deuteron  cross  section  were  measured  under  the  same 
experimental  conditions,  the  deuteron  cross  section  is 

/ 


do 


dfi 


=  N 


dg 
dfi 


fexp 


This  implies  that 


da 


an 


daJ 

dft 


daJ 


deVries 


an 


\ 


da 

dft 


'exp 


/exp 

Now   putting   this    in   terms    of    the    experimental    ratio   R, 
where   R   is    given   by 


R   = 


one   obtains 


da 


an 


'exp 


daJ 


dfi 


'exp 
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da 


dft 


R 


dft 


ydeVries 


and  recalling  that 


da 


dfi 


da  I      r2 

dft    *.+.  d  ' 
'mott 


one  gets  the  final  expression  for  G,  as 


G2  _ 
Gd  ~ 


daJ 


da 


dft   ,  „  .   /  Idfl 

deVries  /  \    /mott 


R   , 
m 


where 


daJ 


dft 


and 


»^>  d 

da 


dfi 
l    'mott 


was  calculated  from  the  BERROTON  program 
deVries 

was  calculated  from  the  BERDEUTR  program  and 


the  ratio  R   is  the  weighted  mean  ratio  given  by  Equation 

IV-17. 

2 

The  statistical  error  in  G,  was  calculated  using  the 

following  equation; 


dof 


deVries 


I*   d 

da 


'R  ' 
m 


dQ,      mott 


where  S    is  given  by  Equation  IV- l! 
m 
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APPENDIX  B 
RMS  STRUCTURE  RADIUS  OF  DEUTERON 

By  substituting  the  deuteron  charge  distribution  for  the 
neutron  charge  distribution,  it  can  be  shown  that  an  expres- 
sion similar  to  Equation  11-15  for  the  deuteron  is 
applicable. 

?            <r^>   5    <r1>   . 
Gd(q  ]    =  Gd(0) 6   q  +  "lfo  g  +  •'•      (B_1) 

2 
where  r ,  is  the  mean  square  radius  of  the  deuteron  charge 

distribution.   In  principle,  the  deuteron  form  factor  could 

2 

be  fit  to  this  polynomial  form  and  the  slope  of  G , (q  )     ~ 

2  2  <rd> 

versus  q   in  the  limit  as  q   approaches  zero  would  be  -  — 

6 
4 
However,  for  the  deuteron,  the  coefficient  of  the  q   term  is 

4 
not  well  determined  and  the  choice  of  <r.>  greatly  affects 

2  4 

the  value  of  <r  >.   In  order  to  determine  <r,>,.  data  at 

d  d 

higher  momentum  transfer  is  necessary  and  its  determination 
requires  knowledge  of  the  coefficient  of  the  q   term. 

Fortunately,  Schumacher  and  Bethe  [7]  have  suggested  a 

procedure  to  avoid  this  problem.   The  method  consists  of 

2 

using  the  experimental  data  and  the  function  C^  (q  )  associ- 

E 

ated  with  a  deuteron  wave  function  to  calculate  A  in  the 
expansion 

Gp  (q2)  =  A2  +  B4  .  (B-2) 


11 


If  the  deuteron  structure  factor  C^  is  also  expanded  for 

2 

small  q  ,  the  expansion  is: 

2         12    2     14    4 
CE«J  »  "  1    -    I  <rE>q   +  120  <rE*  q   +  •••'   (B"3) 

2 

where  <r  >  may  be  defined  as  the  mean  square  radius  of  the 

deuteron  structure  factor. 

Schumacher  and  Bethe  observed  that  if  the  same  data  set 

is  used  to  calculate  A  from  various  C„  functions,  which  have 

hi 

2 

different  values  of  <r  >,  the  plot  of  A  as  a  function  of 

2 

<r  >  forms  a  straight  line,  if  one  uses  only  C„  functions 
hi  hi 

corresponding  to  deuteron  models  which  have  the  correct 
binding  energy.   Then,  if  one  requires  that  A  have  the  value 

deduced  from  the  neutron-electron  scattering  experiment,  the 

2 
value  of  <r  >  may  be  determined,  as  shown  in  Figure  V-8. 
hi 

2 
The  error  m  <r  >  is  fixed  by  the  intersection  of  the  pro- 

h 

jection  of  the  errors  of  A,  and  the  slope  obtained  using 

2 

models  considered  in  this  experiment.   The  value  of  <r  > 

hi 

is  given  in  Section  V. 

Finally,  one  may  note  that  the  mean  square  charge  radius, 

2 

<r ,>  may  be  found  from: 

2       2       2 
<rp    =    <r*>    +    <r*>    _  6A  . 

d       E       p 
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APPENDIX  C 

q4  DEPENDENCE  OF  G 

n 

In  Chapter  IV,  the  slope  of  G    was  determined  assuming 

2      n 
a  linear  fit  of  G    versus  q  ,  but  there  was  some  evidence 

J— i 

n 
that  the  inclusion  of  a  quadratic  term  may  be  more  appro- 
priate.  Therefore,  a  fit  with  both  linear  and  quadratic 
terms  was  made  to  the  data,  but  it  was  discovered  that 
such  a  fit  produced  a  quadratic  coefficient  of  little 
confidence. 

Since  the  value  of  the  coefficient  for  the  linear  term 
is  not  strongly  dependent  on  the  value  of  the  quadratic 
coefficient,  for  this  data,  a  value  for  the  quadratic 

coefficient  was  assumed  using  the  work  of  Galster  et  al . 

4 
[39].   Galster  et  al.  obtained  the  value  B  =  -0.0036  fm 

for  the  quadratic  coefficient  and  this  value  was  used  in 

application  to  this  work's  data.   A  recalculation  for  the 

quadratic  coefficient  of  the  analytic  expression, 

2      4 
GE   =  Aq   +  Bq   , 

n 

shows  that  the  quadratic  term  contributes  5%  to  G^   at 

hi 

2-2  2-2  n 

q   =  0.35  fm    and  10%  at  q   =  0 . 5  f m   .   With  the  assumed 

quadratic  coefficient,  the  slope  of  G^   increased  by  7% 

n 
from  that  obtained  without  the  inclusion  of  the  quadratic 

term.   Also,  the  deuteron  structure  radius  is  decreased  by 

0.15%.   Note  that  the  changes  in  the  slope  and  the  deuteron 
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structure  radius  are  within  the  statistical  accuracy  of 
the  value  previously  reported  in  Chapter  IV.   A  comparative 
table  of  values  utilized  in  the  linear  and  quadratic  fits 
can  be  found  in  Tables  C-l  and  C-2. 
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TABLE  C-l 


COMPARISON  OF  GE   AND  GE  -  Bq 

n        n 


2 

q 

LF5 

LF15 

(fm"2) 

GE 
n 

4 

GE  -B<3 
n 

GE 

n 

GE  -Bq4 
n 

0.05 

0.0002 

0.0002 

0.0000 

0.0000 

0.10 

0.0037 

0.0037 

0.0033 

0.0033 

0.20 

0.0028 

0.0029 

0.0020 

0.0021 

0,25 

0.0070 

0.0072 

0.0061 

0.0063 

0.  30 

0.0076 

0.0079 

0.0066 

0.0069 

0.35 

0.0094 

0.0098 

0.0082 

0.0086 

0.40 

0.0084 

0.0089 

0.0070 

0.0076 

0.477 

0.0099 

0.0107 

0.0083 

0.0091 

0.50 

0.0107 

0.0116 

0.0090 

0.0098 

For  LF1  assuming  B  =  -  0.0036  fm 

A  =  0.0232  ±  0.0022 

4 
For  LF5  assuming  B  =  -  0.0036  fm 

A  =  0.0232  ±  0.0022 

For  LF15  assuming  B  =  -  0.0036  fm4 

A  =  0.0198  ±  0.0023 


TABLE  C-2 

COMPARISON  OF  THE  DEUTERON  STRUCTURE  RADIUS  OBTAINED 
FROM  THE  LINEAR  AND  QUADRATIC  FITS  TO 
G^  ,  ASSUMING  GROSS  RELATIVISTIC  CORRECTION 
n 


G    =  Aq' 


2     4 
G    =  AqZ  +  bq 

En 


<4> 


3.867  ±  0.018 


3.858  ±  0.011 


rms 


1.9665  ±  0.0045 


1.9642  ±  0.0045 


APPENDIX  D 
COMPARISON  OF  GROSS1  AND  FRIAR'S  RELATIVISTIC  CORRECTION 

Both  the  Gross  and  Friar  relativistic  corrections 
multiply  the  right  hand  side  of  Equation  II-7.   The  Gross 
relativistic  correction  factor  is  given  by: 


_  J3L 


=  1  - 


2 
-9-  + 


8M 


2  '      4  ' 
4M     64  M 
P        P 


(D-l) 


and  Friar's  correction  factor  is 


4M' 


Expanding  Friar's  correction  in  a  binomial  series  to  second 

2 
order  in  3_.  f    the  following  is  obtained: 


4M 


1  + 


-1 


SL 


4M' 


2       4 

=  1  -  L  +  _SL_ 


4M' 


16M 


4  ' 


(D-2) 


The  difference  between  the  two  relativistic  correction 
factors  is  approximately  given  by  the  difference  between 
Equation  D-2  and  D-l.   If  this  difference  is  given  the 
symbol  A  then, 


A  -  SL 


M 


_1_ 
16 


_1_ 

64 


=  2.    2l_ 

64    M   4   ' 

M 


(D-3) 


The  difference  is  seen  to  increase  with  increasing  q  . 

2-2  2 

At  q   =  0.5  fm  ,    which  is  the  highest  q   investigated,  A 
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-5 

is  2.3  x  10    with  Friar's  correction  always  being  greater. 

2         -2 
Since  the  corrections  are  both  about  0.997  for  q  =  0 . 5  f m   , 

then  the  difference  in  these  two  corrections  is  about  two 

5 
parts  in  10  .   Table  D-l  gives  a  comparison  of  both  correc- 

2 
tions  for  all  q  's  investigated 

Friar  [8]  also  shows  that  C^ ,  the  electric  structure 

Ei 

2 
factor  of  the  deuteron,  should  not  be  evaluated  at  q   but 

2' 
at  an  effective  q    given  by 

q    =  SL- .  (D-4) 

i  +  a 


4M^ 
d 


This  also  creates  a  very  small  change  as  seen  in  Tables  D-2, 
D-3  and  D-4.   Friar's  first  relativistic  correction  opposes 
the  tendency  of  the  use  of  the  effective  momentum  transfer, 

so  the  net  difference  between  the  use  of  the  Friar  and  Gross 

4 
corrections  is  smaller  than  four  parts  in  10  .   Also 

included  in  these  tables  is  a  comparison  of  the  slope  of 

2 

G   versus  q   using  each  correction.   Again  the  difference 

n 
is  seen  to  be  negligible  for  this  experiment. 
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TABLE  D-l 
COMPARISON  OF  RELATIVISTIC  CORRECTION  FACTORS 


q2 
fm"2 

Friar ' s 

Gross ' 

,   Friar 
Gross 

Relativistic 
Correction 

Relativistic 
Correction 

0.05 

.999972 

.999972 

-.000000 

0.10 

.999447 

.999448 

-.000001 

0.20 

.998894 

.998896 

-.000002 

0.25 

.998618 

.998621 

-.000003 

0.30 

.998341 

.998345 

-.000004 

0.35 

.998065 

.998070 

-.000006 

0.40 

.997788 

.997796 

-.000007 

0.477 

.997360 

.997370 

-.000010 

0.50 

.997235 

.997247 

-.000011 

TABLE  D-2 

2  2' 

COMPARISON  OF  C„  (q  )  AND  C„  (q   )  FOR  LFl 


2 

q 

CE(q2) 

2 

q 

CE(q2') 

CE(q^  ) 
CE(q2) 

0.05 

.96882 

0.0500 

.96882 

-.00000 

0.10 

.93967 

0.1000 

.93969 

-.00001 

0.20 

.88668 

0.1999 

.88673 

-.00006 

0.25 

.86247 

0.2498 

.86256 

-.00009 

0.30 

.83961 

0.2998 

.83972 

-.00013 

0.35 

.81796 

0.3497 

.81810 

-.00017 

0.40 

.79742 

0.3996 

.79760 

-.00023 

0.477 

.76766 

0.4765 

.76790 

-.00030 

0.50 

.75929 

0.4997 

.75955 

-.00034 

Slope  of  Gj?   versus  q   using  Gross'  correction  is 

n      0.0218  ±  0.0022. 
Slope  of  GE   versus  q2  using  Friar's  correction  is 


■n 


0.0213  ±  0.0022. 
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TABLE  D-3 

2  2' 

COMPARISON  OF  C„ (q  )  and  C^ (q   )  FOR  LF5 

hi  hi 


2 

q 

CE(q2) 

2 

q 

CE(q2') 

,   CE«*2'> 

CE(q2) 

0.05 

.96881 

0.0500 

.96881 

-.00000 

0.10 

.93966 

0.100 

.93968 

-.00002 

0.20 

.88665 

0.1999 

.88670 

-.00006 

0.25 

.86244 

0.2498 

.86252 

-.00010 

0.30 

.83956 

0.2998 

.83968 

-.00014 

0.35 

.81790 

0.3497 

.81805 

-.00018 

0.40 

.79735 

0.3996 

.79754 

-.00023 

0.477 

.76759 

0.4765 

.76782 

-.00031 

0.50 

.75921 

0.4997 

.75947 

-.00034 

Slope  of 

GEn 

2 

versus  q^ 

0.0220 

using  Gross 
±  0.0022. 

'  correction  is 

Slope  of 

GEn 

2 
versus  q 

0.0214 

using  Friar 
±  0.0022. 

' s  correction  is 

TABLE  D-- 

4 

COMPARISON  OF  C^(q2)  AND  C„(q2  )  FOR  LF15 


q2 

CE(q2) 

2 

q 

CE(q2') 

.   CE(02'> 

2 
CE(q2) 

0.05 

.96900 

0.0500 

.96900 

-.00000 

0.10 

.94002 

0.1000 

.94003 

-.00002 

0.20 

.88731 

0.1999 

.88737 

-.00006 

0.25 

.86324 

0.2498 

.86332 

-.00010 

0.30 

.84049 

0.2998 

.84060 

-.00013 

0.35 

.81895 

0.3497 

.81909 

-.00018 

0.40 

.79851 

0.3996 

.79869 

-.00023 

0.477 

.76889 

0.4765 

.76913 

-.00030 

0.50 

.76056 

0.4997 

.76082 

-.00034 

Slope  of  Ge versus  q^  using  Gross' correction  is 

n      0.0185  ±  0.0022. 
Slope  of  Ge   versus  q2  using  Friar's  correction  is 

n      0.0180  ±  0.0022. 
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APPENDIX  E 
COMPUTER  PROGRAMS 

Six  FORTRAN  programs  were  used  to  perform  calculations 
needed  for  this  work.   They  are  TOPRADCR,  TOPROTON,  TOPDEUTR, 
BERHFORM,  TRADFORM  and  BERRELCR.   The  first  three  were 
written  by  Topping  and  can  be  found  in  Reference  [5] . 
BERHFORM  and  TRADFORM  calculate  the  experimental  cross 
sections  for  each  of  the  four  channels  considered  for 
hydrogen  and  deuterium  respectively.   BERHFORM  is  repro- 
duced on  the  following  pages.   BERRELCR  was  used  to 
compare  the  relativistic  corrections  of  Gross  and  Friar. 
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